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ABSTRACT 


Stereology is a collection of mathematical methods that relate the parameters 
that define a spatial structure to lower dimensional measurements made on sections 
through the structure. Applications of such methods are found in various sciences, such 
as materials science, medicine, biology, forestry, geology etc. Quantitative 
characterization of materials microstructures often require estimation of hue 
particle/grain size distributions. In this study, numerical methods are successfully 
developed for de-convolution of apparent 2D profiles of spheres, variable spheroids and 
circular disks. It is found that use of profiles of spheroids as equivalent area circles to 
de-convolute 3D size distribution is not a bad approximation. Effect of sample size, 
width of class interval and mode of representation of mean size of a bin, on the 
estimates of distributions and mean particle sizes is studied. It is realized that using root 
mean square value (in de-convolution of spheres) to represent mean bin size provides 
the most bias firee estimate of mean particle size. 
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CHAPTER 1 


INTRODUCTION 

Stereology (also known as quantitative microscopy or quantitative metallography) 
is the area of stochastics dealing with statistical inference about spatial structures from 
geometric samples of the structure such as two-dimensional sections and one- 
dimensional probes. It can be described as geometric sampling theory or a set of 
sampling methods for quantifying geometric parameters such as the length of 
dislocations, volume fraction of a particular phase, or average size of grains, in a given 
microstructure, or the volume of an organ. Stereology provides methods for sampling 
from such structures and thereby for estimating the parameter of interest. The 
parameter can often be estimated from simple measurements of the samples, e.g. 
counting. The development of stereological methods involve the use of advanced 
mathematical tools, especially from geometric probability theory and integral geometry 
[1-2]. Stereology may be applied in various sciences, such as materials science, 
medicine, biology, forestry and agriculture. The word 'stereology" originates from the 
Greek word 'stereos': solid, - here in the sense of having all the geometrical dimensions 
or being cubic [3]. So stereology traditionally concerns three-dimensional structures. 
However stereology is also applicable to structures that are subsets of a two- 
dimensional space. 

Physically cutting a specimen of material and looking at the distribution of the 
particle intersections on a lower dimensional section usually obtain the raw data. The 
estimation of properties of the particulate phase (‘particulate phase’ is meant to include 
voids i.e. lack of material as well as true particles) using such raw data splits into two 
broad categories [4]. The first is the estimation of aggregate properties of the particulate 
phase, which are definable as average aggregate properties of particulate phase 
intersections on a lower dimensional section. Probably the most familiar examples of 
this type of analysis are area and lineal analysis volume fraction estimates. 

The second category of estimation procedures depends upon first obtaining an 
estimate of the particle size frequency function. Moments of the estimated frequency 
function are then used to estimate functions of moments, such as mean particle size and 
surface to volume ratio. These estimation problems are not peculiar to metallurgy. 
Clearly, the particulate phase could be nebulae in outer space, leukocytes in blood. 
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aggregate in cement, nucleoli in nerve ganglion or trace mineral in rock as opposed to a 
metallurgy situation, such as carbon particles dispersed in steel. 

In general the estimate depends upon particle shape and orientation. A system of 
particles is often represented simply by the mean particle diameter or by the particles 
per unit volume. However in order to characterize a particulate structure more 
completely, knowledge of the distribution of particle sizes is required. Microstructures 
of materials often contain particles/grains/voids/inclusions having a range of sizes and 
shapes. 3D size and shape distribution is an important attribute of such microstructures, 
since the mechanical behavior of these materials largely depend on characteristics of 
particle/grain size distribution. In such processes as crystallization from the melt, 
recrystallization and nucleation and growth of precipitate particles, it may be desired to 
know not only the size distribution, but also how this distribution changes with time. If 
the microstmcture consists of population of only spherical particles having different 
size, then it can be characterized by a simple univariate size distribution function. But 
microstructures often contain a population of particles having different sizes and shapes 
and in such cases, more than one size parameter is required to uniquely identify any 
given particle in the population. Such microstructures are characterized by bivariate or 
multivariate size distribution functions. Global parameters such as the number density 
(Nv) of particles, are not sufficient for the microstmcture characterization during 
different processes of grain formation and grain growth. It is necessary to describe, the 
population of particles by parameters corresponding to their individual geometry, 
namely the size, shape factor and orientation. However, these may interact with each 
other and, therefore, the joint probability distribution of geometrical parameters has to 
be estimated. 

Present study aims at developing numerical procedures for the estimation of 3D size 

distribution of spherical particles, bivariate size-shape distribution of variable spheroids 

and bivariate size-orientation distribution of circular disks from the corresponding 2D 

profile measurements. The appropriateness of assuming spherical shapes for non- 

spherical systems to de-convolute 3D distributions is studied. Effect of variables like 

sample size, class interval and . choice of mean bin size on 3D estimations is also 

analyzed. In addition to this introductory chapter this thesis consists of five chapters. In 

second chapter a detailed review of the development of different stereological methods 

over the past years is given. Whereas, third chapter consists of the theoretical 

background of the stereological methods for unfolding the probability distribution of 
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geometrical parameters of spheres and spheroids. The results obtained in the present 
study are given in fourth chapter along with a brief discussion about the results. The 
fifth chapter consists of the concluding remarks, followed by a chapter on suggestions 
for future work. 
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CHAPTER 2 


BACKGROUND 

As mentioned in the introduction, stereological estimation methods fall into two 
broad categories: 1) Methods for estimation of aggregate properties of the particles, 2) 
Methods for estimation of distribution function of geometric parameters of the 
particles. A brief overview is given in this section about the development of these 
methods and their underlying principles. 

2.1 METHODS FOR ESTIMATION OF AGGREGATE PROPERTIES 


The earliest work on quantitative microscopic analysis goes back to French 
petrographer A. Delesse (1848) , who proved mathematically that in a fully random 
cross section of a uniform aggregate, the area occupied by each constituent of the 
aggregate is exactly proportional to its volume in the mass of the rock [5]. Delesse 
measured the relative areas of minerals lying in a macroscopic section by cutting out 
segments of thin foil to represent the areas occupied by the various minerals. He then 
weighed the pieces of tin corresponding to a given mineral, and divided this by the total 
weight of tin corresponding to the entire structure. Sorby (1856) extended the areal 
method of analysis to microstructures by using a camera for tracing on “an evenly thick 
piece of drawing paper” the outlines of minerals seen imder a petrographic microscope 
[6]. He weighed tbe whole sketch, and then cut the areas representing the various 
minerals and weighed them as Delesse did with tin foil for much coarser structures. 
Sollas (1887) combined the two methods by cutting out areas of tin foil corresponding 
to a camera tracing [7]. Julien (1903) made similar determinations in a petrographic 
study, but attained greater precision by photographing the structure, and then cutting 
and weighing the paper segments of the micrographic print [8]. A few years later, 
Rosiwal (1898) showed that the volume proportions of minerals in uniform rock are 
equal to the lineal proportions intercepted by the minerals on a random line passed 
through the structure [9]. He not only demonstrated the applicability of lineal analysis 
to geometrically produced areas, but also found good agreement with the areal method 
of Delesse as applied to actual macroscopic sections. For microstudies, Rosiwal 

employed a micrometer ocular to measure the fraction of a line-traverse intersected by 
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the minerals of each type. A number of areal and lineal methods of measurement have 
been developed with a view to increasing precision and decreasing tedium. Joly (1903) 
traced the enlarged image of a microstructure (projected onto the ground glass of a 
camera) on graph paper, and estimated the proportions of the individual squares 
occupied by the various minerals [10]. He then summed up the total area of each 
microconstituent in the polished section. Johaimsen (1919) made tracings of 
microconstituents with the aid of a camera, and measured the areas with a regular 
planimeter [11]. The earliest applications of point ratios as estimators of areal ratios 
appear to have been made by Thomson (1930). He used one and two-dimensional point 
grids in order to determine the fractional number of points falling within the boundaries 
of a constituent in the section plane [12]. Lightner and Herty (1932) improved upon the 
point-counting procedure by photographically superimposing a grid on a print of the 
microsfructure [13]. Gardner (1943) carried out point-counting by covering the 
micrographic prints with translucent graph paper [14]. H. S. Rawdon proposed a 
photometric method for determining relative areas of dark and light constituents in a 
two-constituent structure but it is hardly suitable for routine work [15]. Polushkin 
(1925) appears to have been the first to apply lineal analysis to metallic structures (Bi- 
Sn eutectic) [16]. He mounted micrographs on a movable platform, and used a 
recording tape to sum up the lengths of a given constituent intercepted during a 
traverse. He was able to determine the proportion of eutectics and eutectoids in various 
alloys by this method of analysis. A number of critical petrographic studies have been 
made to compare the areal and lineal methods of quantitative analysis. Thomson (1930) 
has applied the Delesse, Rosiwal and grid methods to synthetic aggregates and actual 
stmctures, and has found agreement of ail three techniques within a few percent. 
Karlsson and Drar (1980) demonstrated the volume fraction determination by the aid of 
computer modeling [17]. In the nomenclature of quantitative metallography developed 
by Underwood [18] this is expressed by 
Equation 2-1 

Vy=A,=L,=P, 

Saltykov (1945) gave another major breakthrough m stereological estimation 
methods, when he first derived the following basic equation for obtaining the area of 
surfaces in a volume, just from a simple counting on two-dimensional section [19]. 
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Equation 2-2 


Sy = 2P^ 


The importance of this equation can be appreciated from the fact that it has been 
rederived at least eight times since 1945. This equation applies to any system of 
surfaces in space, whether plane or curved, continuous or interrupted, isolated or 
connected i.e. systems of surfaces with any configuration. It is as valid for systems of 
interconnected surfaces as for systems of discontinuous, separated, or bounded 
surfaces. Tomkeieff (1945) proposed a relationship between the surface-to- volume ratio 
of an isolated, convex body and the mean intercept length through the body [20]. 
Equation 2-3 



L, 


It was later claimed, however, that the requirement of convexity is unduly restrictive. 
Chalkley (1949) devised a method for obtaining the average surface-to-volume ratio of 
a system of particles, which involves both a point count and an intersection count on 
random section through the particles [21]. 

Equation 2-4 



Assuming that all the grains are the same size, Kaiser expressed Sy , the intergranular 


area per unit volume, in terms of Ny , the number of grains per unit volume [22], by 


Equation 2-5 


Sy = FN^ 

Where F is a constant dependent on geometric shape. Values of F are 3 for cube, 2,866 
for the hexagonal prism, 2.673 for the rhombic dodecahedron, and 2.659 for the 
truncated octahedron. For real grains, which were thought to approach the properties of 
rhombic dodecahedra most closely, Kaiser’s equation is approximated as 
Equation 2-6 





Guliver (1918) derived the following simple expression relating the number of grains in 
space and the planar grain area [23] 
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Equation 2-7 



Combining the above two equations gives the following expression 
Equation 2-8 


I 



In an attempt to approximate the real structural shapes of polycrystalline aggregates 
more closely, some investigators abandoned the assumption of equal size and identical 
shape of all the constituent crystallites. Meijering (1939) determined analytically the 
parameters of the spatial structure of a metal polycrystal in which the crystal start to 
grow simultaneously and isotropically from nuclei distributed at random [24]. For the 
specific area, he derived the following relationship as a function of the number of 
grains per unit volume: 

Equation 2-9 

1 

Sy = 2.9 

Approximate methods have been proposed for obtaining the specific surface area of 
discrete particles, provided that one can assume they have diameters of the same size. If 
the particles are spheres of diameter D, and the number of particles per unit volume is 
Ny, then their specific surface will be 
Equation 2-10 

Sy =7ZD^Ny 

Volume fraction of the particles will be 
Equation 2-11 

_nD^Ny 

6 

From the above two equations eliminating Ny yields 
Equation 2-12 
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S. Z. Bokshteyn developed the above equation by a different method, and then replaced 

D by the mean diameter of the particle cross sections d [25]. This gives the specific 
surface area of a monodispersed system. 

Equation 2-13 

d 

with an error not more than 7 percent. Instead of making the approximate substitution. 


d zD, one may use Pullman’s statistically exact relation for uniformly sized spheres 
[26]. 

Equation 2-14 



This substitution yields 
Equation 2-15 



In order to account for a distribution of sizes in the 3-D particles, Saltykov (1954) 
applied a correction term to the above relation. 

Equation 2-16 




6Vy f TT 1 ^ 


Where S is the ratio of <y{D], the root-mean-square deviation of the diameter of the 


spatial particles to D , their mean diameter, the coefficient of variance of the particle- 
size-distribution function. In order to calculate the surface area of spheroidal carbides, 
L. S. Moroz (1950) developed the expression [27] 

Equation 2-17 

Sy =1.68VQW^ 

Where Ck is the ratio of the volume firaction of the carbides to their weight percent of 
carbon (6.67). Cahn and Nutting (1959) investigated particle characteristics in thin 
films by electron transmission microscopy [28]. They derived the following expression 
for the surface area of these particles 
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Equation 2-18 


Sy =^[^A{proj)-Vy] 

Where Aa (prof) is the projected area fraction of convex particles, and t is the thickness 
of the thin film. Their derivation utilized the following general expression for the 
surface area S of convex bodies, attributed to Cauchy [29] : 

Equation 2-19 

S = AAiproj) 

A familiar constituent of alloys is the lamellar structure. Perhaps the best-known 
example is pearlite in steels. The interface area of a lamellar structure depends on the 
magnitude of the spacing between lamellae. Gensamer, et al. (1942) proposed a method 


to calculate the mean intercept spacing 7 [30]. Many straight secants of definite length 
are drawn randomly across the microstructure, intersecting a great number of platelets 
at all possible angles. The total length of test secants divided by the total number of 
cementite platelets intersected is defined as the mean intercept spacing. Combining this 
definition and the basic equation for specific surface area they got the following 
expression 
Equation 2-20 



In 1953, C. S. Smith and L. Guttman developed a new stereological method to 
find the length of lineal elements on a plane from a simple counting technique and 
derived the following expression [31]. 

Equation 2-21 



They also framed a stereological technique to find the total length of random 
lines in a given volume from a simple two-dimensional measurement. Mathematically 
their technique can be expressed as. 

Equation 2-22 

Ly =1P^ 
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An equation for the ratio of perimeter length to area of isolated two-dimensional 
convex figures was given by Tomkeieff (1945), which he expressed in terms of the 
mean intercept length. 

Equation 2-23 



A £' 


It would be possible to observe true points in a volume only under unusual 
circumstances. The intersections of dislocations in planar arrays are observed 
frequently in thin metallic foils by means of electron transmission. The intersections 
actually represent a set of points in a three dimensional volume. Williams and Smith 
(1952) from their experiments established that if the grain edges of a polycrystalline 
aggregate were delineated, one could observe a three-dimensional network of lines 
(edges) terminating four at a grain comer (point) [32]. Milliards (1962) developed a 
method to obtain no of points per unit volume in space [33]. In mathematical terms his 
method can be expressed as, 

Equation 2-24 

Py = ^PaPl 

Gurland (1958) derived an equation to find the contiguity C, which he defined 
as the firaction of the total interface area of a phase that is shared by particles of the 
same phase [34]. For a particulate stmcture with particles of a-phase in a homogeneous 
matrix of p-phase, he defined the contiguity as 
Equation 2-25 

^ <PlL 

Where (F,L and (F,L are the number of intersections per unit length of test line 

with aa and ap interface traces. Cahn and Hilliard (1959) proposed other parameters 
for describing quantitative aspects of particle contact [35]. One such parameter is 
(Ly , the aap perimeter length per unit volume around contact areas between a-a 

particles; another parameter is {ly , the aaa edge length per unit volume between 
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three contiguous a-particles. The total perimeter length per unit volume around aa- 
surfaces, (Ly , consists of contributions from both aaa- and aaP-lines. Hence, the 

following equation is apparent from the above discussion: 

Equation 2-26 

(^yL‘6(pJ^+2(pJ^ 

Where the coefficient 6 includes a factor of three because each aaa-edge is shared by 
three aa-contacts. Each aaP-edge, on the other hand, adjoins only one aa-contact. Cahn 
and Hilliard (1959) derived an expression for (■Ly)^^ in terms of the numbers of 

interface traces seen in a unit area on the section plane. Based on the above mentioned 
fact they obtained the following equation from which the final equation is obvious: 
Equation 2-27 




'aa/S 


Thus the perimeter lengths (ly can be obtained either by counting contact traces or 

by counting “triple” points. The factor H in the first of the above two equations arises 
because each interface trace has two termini on the section plane. The above 
expressions dealing with contiguity are perfectly general, and can be evaluated by 
simple measurements on the section plane. Gurland (1958) has defined another quantity 
average number of contacts per grain (N^. by incorporating the shape factor and the 


average and standard deviation of perimeter lengths as below: 
Equation 2-28 


(Nc). 



Where the term in the numerator represents the number of a-grain contacts and the 
denominator is the total number of a-grains, both per unit volume. It is apparent that 
this expression can’t be evaluated by measurements on a single section plane unless a 
particle shape is assumed. Forscher (1955) defined another attribute of particulate 
system, called the “continuity” of the particles, as the probability of the presence of 
infinitely long chains of connected particles [36]. Consideration of the probabilities of 
such chain formation fall under the theory of branching chains that is too far afield. 
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In 1967 J. W. Cahn and R. T. Dehoff independently developed a method for 
estimating the mean mean curvature of a spatial surface by a tangent count on plane 
sections [37-38]. Their formula is: 

Equation 2-29 

k = -C = -TJL 
4 2 ^ ^ 

where is the number of intersections per unit test line length, and is the net 
tangent count per unit section area. 

Oriented structures occur so frequently in metallic alloys that special procedures 
have been developed to determine their properties. Meteorites, twinned minerals or 
metals, banded rocks, dislocation arrays, nonmetallic “stringers” in extruded rods and 
the lamellar eutectic alloys are a few of the examples. The Stereological equations 
mentioned above are derived for random structures. For oriented systems analogous 
relations hold, but with changed coefficients. 

The degree of orientation of a system of lines in a plane may be quantitatively 
defined as the ratio of the length of the oriented portions of the lines to the total length 
[39]. To estimate this property of a structure, it is necessary to determine separately the 
oriented portions of interest and compare their length to the total length of lines in the 
system. The oriented portion of the lines can be calculated separately from the 
isometric portion by measuring Ni at certain directions relative to the orientation axes. 
In a system with one orientation axis, two groups of secants would be used- one parallel 
to the orientation axis, the other perpendicular to it. This method has been called 
“directed secants on a plane” and is proposed by Saltykov. A partially oriented system 
of lines may be regarded as consisting of two superimposed systems. Then the 
numerical expression for the degree of orientation ©of a system of lines will be 
determined, in percent, by: 

Equation 2-30 

where is the specific length oforiented portion of lines 

(LyOij is the specific length of isometric portion of lines 
Saltykov proposed a special polar plot to summarize graphically the orientation of all 
line segments in the system of lines. Test lines are uniformly distributed over all 
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possible angels to the orientation axis, and then PL(d) is obtained for each angle. By 
using polar coordinates, the angles of the lines and the specific number of intersection 
can be represented readily by a curve, which he called as the “rose of the number of 
intersections” or simply “rose.” Underwood derived a relation to find the length per 
unit area of a system of oriented (parallel) lines. 

Equation 2-3 1 

But in real stmctures the lines in a microsection are usually either isometric or partially 
oriented. Only in rare cases, such as the Widmanstatten patterns in meteorites, do we 
find a completely oriented system of lines. In a partially oriented system of lines in a 
plane, part of the total length of lines is oriented in a definite direction (or directions). 
The remaining segments may have essentially a random orientation. Underwood 
derived another equation to find the La in such cases. 

Equation 2-32 

(lX =1 (pJ, 

Now the degree of orientation of the lines in the system can be expressed as, 

Equation 2-33 

Completely oriented systems of lines in space with one orientation axis consist of 
straight lines that have equal or variable spacing, but are parallel to each other and to 
the orientation axis. Underwood derived an equation to find the length per urdt volmne 
of such oriented system of lines in space. 

Equation 2-34 

The coefficient in the above equation is unity, whereas when the test planes are applied 
randomly to a system of lines in space, the coefficient is two. Schoeck (1962) 
investigated various planar arrays of dislocation lines in space and obtained coefficients 
with values between one and two, depending on the relative orientation of the section 
plane with respect to the dislocation arrays [40]. In case of partially oriented systems of 
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lines in space an equation analogous to the one described above for a partially oriented 
system of lines in plane is derived. 

Equation 2-35 

He has derived the following equation to find the siurface area per unit volume of the 
oriented system of plane surfaces. 

Equation 2-36 

\r ~ )ii 

In case of partially oriented system of surfaces he obtained the following equation. 
Equation 2-37 

5,=1.57l(P,)^+0.429(PJ, 

The total surface Sy is equal, in general, to 
Equation 2-38 

Sy = {Sy + (Sy + {S y ) 

Where, is, lin, pi refer to isometric, linear, and planar components of the systems of 
surfaces, respectively. The general equation for the degree of orientation of planes in 
space, in percent is 
Equation 2-39 

lOOfo)., 

^or ^ 

Milliards (1967) provided a new mathematical firamework for the quantitative 
description of anisotropic microstructures [41]. He has developed a self-consistent and 
generalized mathematical treatment for both specifying and determining anisotropy. 
His approach is a completely general one. Two systems of lines are superimposed on 
one another and their probabilities of intersection are considered. One system of lines is 
arbitrarily designated as the test array, and the general equation is set up relating the 
line length of the other system of lines to the ntimber of intersections with the test 
array. Methods are then indicated for obtaining the specific line length as a function of 
direction. Similar procedures are followed for systems of surfaces. The advantages of 
his analysis are that the functions are easy to interpret and provide for direct correlation 
of structural anisotropy with other properties; disadvantages lie in the number of 
angular measurement required and in the mathematical complexity of the solutions. 
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Spektor (1954) has also given an exact treatment for oriented systems of surfaces, but 
his method is restricted to symmetrical structures [42]. The equations derived by 
Hilliard for Sy, La, LvZxq given below. 

Equation 2-40 

La 

Ly = lKLy{^,d) 

Sy =27tSy{{f/,a}) 

p = — Sx,L„ 
y 2 

Pullman (1953) derived the following stereological relation to obtain a spatial 
parameter of great importance in particulate systems, the mean free distance, L, 
between particles, say, of a-phase [26]. 

Equation 2-41 

This equation is valid regardless of size, shape, or distribution of particles or other 
separated volume elements. Chandrasekhar (1943), Gurland (1961), Ashby (1966), 
Westmacott (1966) [43-45] defined some other distance and spacing parameters for 
expressing separation of particles. Depending on the desired application of the 
particular “spacing” the definitions are based on one, two, or more “nearest” -neighbor 
particles. Gurland used two precisely defined nearest-neighbor distances. These two 
distances are denoted by A 2 and A 3 and are defined as the average distance between 
nearest-neighbor pairs of point particles in a plane and in a volume, respectively. 
Equation 2-42 

A. =0.5P/>^ 

A 3 =0.554P/>^ 

Frequently the properties of a particulate or cellular structure can be represented 
satisfactorily by calculations based on a regular polyhedron or voliune of revolution. 
Coexter (1963), Lyustemik (1963) have contributed to the geometrical analysis of 
regular polyhedrons [46-47]. Tomkeieff (1945), Hadwiger (1955), Kendall and Moran 
(1963) have dealt with general relationships for convex particles and figures of 
unspecified shape that are basic to a geometrical analysis. Smith and Guttmann (1953) 
presented general relationships for two-dimensional figures and three-dimensional 

objects, which assumes significance in the present context. Rutherford, Abom, and 

15 



Bain (1937) made calculations on the tnuicated octahedron that form the basis for the 
ASM Metals Handbook tabulation of spatial and planar grain properties [48]. Kaiser 
(1938) tabulated properties of four basic space-filling regular polyhedrons with plane 
faces: i.e., the cube, hexagonal prism, rhombic dodecahedron, and the truncated 
octahedron. Hull and Houdk (1953) selected four regular polyhedrons for an analysis of 
intercept area distribution curves [49]. Saltykov (1958) reviewed his work in this field, 
and contributed some interesting analyses for spheres. The sphere has been substituted 
for particles or cells more than any other volume of simple shape. If it can properly be 
assumed that spheres approximate the real stmcture satisfactorily, one can readily 
obtain spatial parameters firom simple planar measurements alone. The problem is still 
simpler if a constant sphere diameter can be assumed. Pullman (1953) showed the 
power of statistical-geometrical probability approach to the analysis of simple particles. 
He developed methods for the analysis of particle size measurements by area and lineal 
analysis when the particles are in the form of spheres, plates, or rods for both 
monodispersed and polydispersed systems of particles [26]. He derived quantitative 
relationships between various average quantities obtained from a polished cross section 
and the actual dimensions of the particles present. In case of spheres he has derived the 
following relations based on the fact that the probability of a plane intersecting a single 
sphere positioned at random within the unit cube is 2r. This fact is apparent since, of 
the various possible positions for the cross-sectional plane over the unit length from top 
to bottom of the cube, only those positions existing over the length 2r would lead to the 
plane intersecting the sphere. 

Equation 2-43 

A-=—nr 

3 

Nj^ =Ny7tr^ 

1 ^ 

/ =—r 

3 

He derived similar relations for particles in the form of circular plates of iiniform radius 
(radius r, thickness t, r»t). He proved mathematically that the probability of a plate 
intersecting a singly randomly positioned and randomly oriented plate is equal to idl 
times the radius of the plate. Based on this relation he derived the following equations: 
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Equation 2-44 


A = 2rt 


l = 2t 

He also dealt with the particles in the form of randomly oriented rods (radius r, length 
H, H»r). He showed that the probability of a single randomly positioned and oriented 
rod cut by a plane is equal to half of the height of the rod. Based on this probability he 
derived the following equations for rod shaped particles: 

Equation 2-45 


A = 27 cr~ 






nrH 

2 


I = 2r 


It can be seen here that neither the average cross-sectional area nor the average lineal 
traverse of particles in the form of long thin rods depends on the length of the rods. 
This is due to the fact that for such particles there is negligible probability that a cross- 
sectioning plane will cut along the rod length. Later on his analysis is extended for 
many other shapes of the particles. Dehoff (1961) has treated ellipsoidal (prolate and 
oblate) particles in general [50], where as Myers (1963) treated simple particles with 
flat sides (regular polygons) [51]. All the equations given above are for monodispersed 
system of particles. Mean values of V,S,Nv and r for polydispersed systems of spheres 
and think circular disks are derived by Pullman. He divided the particles into class 
intervals, each of constant size, applied the equations described above to each class 
interval, then summed up all valued over the range of particle sizes. However, another 
measurement is required in order to express the average values in terms of measurable 
two-dimensional quantities. For a distribution of spheres, the required quantity is the 

mean value m of the reciprocals of the circle diameters on the test section. Some of the 
equations obtained by Pullman for spheres having a distribution of sizes are reproduced 
below: 
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Equation 2-46 


V = rcVyl2mN^ 

S = 27iN^/mN^ 

Ny = 2m / 7t 

r = 7 t ! Am 

Pullman also formulated the equations for polydispersed systems of thin disks. They 
are given below: 

Equation 2-47 

V = 7V^VylZEN^ 

S = 7t^NJ2EN^ 

Ny=%ENj7C^ 
r = 7vlAE 
t = 7tGI AE 

/ 

Where E is the mean value of the reciprocals of the disk lengths on the section plane 

and Gis the mean value of the ratios of the width-to-length dimensions of disk 
sections. 

The number of sections seen on a polished cross-section and the actual number 
of profiles in the volume are related by the following relation, which is treated as one of 
the basic stereological relations. 

Equation 2-48 

N,=NyH' 

Where H is the mean of the projected heights of the system of randomly oriented 

particles. This relation is exact for convex particles of any shape. But H can’t be 
determined directly in many cases. For a monodispersed system of spheres, how ever, it 

is clear that the diameter D is equal to i/' and in a polydispersed system of spheres, D 

is equal to H . Thus, for spheres of equal size, the above equation becomes 
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Equation 2-49 


Ny = Njb 

Weibel and Gomez (1962) found the following relationship between the numerical 
density of particles, Ny , the numerical profile density, , and the volume finction of 

particles, as a measure of particle size [52]; 

Equation 2-50 


Where, the coefficient J3 is solely dependent on the shape of the particles and not on 


their size, and is related to the voliune v and average profile area of the particle a as 
Equation 2-51 

£ 

V = /3.a^ 

Sterio(1984) proposed an unbiased estimator of number and size of aribitraxy 
particles in 3D space, known as disector. By centering on the design of sampling 
strategies, the disector method first dropped all inherent assumptions about the particle 
shape in previous analysis methods, then the point-sampled intercept and some other 
methods followed. The disector comprised, for the first time, a complete rule for 
uniform sampling of 3 -dimensional objects and other discrete events according to 
cardinality [53]. The disector comprises an integral test system and a counting rule: the 
number of structures sampled by the counting firame disappearing firom one section 
plane to another, Q- is counted. The integral test system comprises the 2-dimensional 
counting firame with a certain area per test point, (a/p), and the number of times the test 
points hit the reference volume under study, P. The two parallel section planes are a 
known distance, h, apart and the numerical density of stmctures is estimated as: 

Equation 2-52 

/V = 

'' {a/p)Jt.P 

Sampling with the disector is not restricted to isolated objects. Associate features of 
objects may also be sampled by the disector. Disector sampling has opened up for a 
whole range of local stereological estimators of size. However if the matter is opaque. 
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the physical dissector has to be used which is based on parallel sections a known 
distance apart. Then it is necessary to identify particles that are hit by both sections, 
which requires perfect registration of the relative position of the two images. Moreover 
a bias may be introduced because particles that are smaller than the distance between 
the two section planes cannot be detected. Often it is practically very difficult or even 
impossible to make parallel section, as of hard and brittle material. This is a reason why 
classical stereo logical methods using single sections are traditionally preferred by 
material scientists. 

2 . 2 METHODS FOR ESTIMATION OF DISTRIBUTION 
■ FUNCTIONS 

Estimation of the distribution of grain sizes in three dimensions in space is 
another similar problem that is of interest to metallurgists. Although the earliest work 
on quantitative microscopic analysis goes back to Delesse (1848), the estimation of 
particle size frequency functions dates only from Wicksell (1925) that is in connection 
with the biometric examination of corpuscles [54]. Through an integral equation, he 
related the size frequency function for spherical particles to that of circular intersections 
on a random plane of intersection. The fundamental treatise of Wicksell is one of the 
masterpieces of stereology. 

Many different approaches have been used to solve Wicksell's problem or 
variants of it. The oldest technique for finding a solution is by discretizing the 
distribution of the sphere radii and the distribution of the planar circular profiles, the 
last discretization being estimated by a histogram of the observed data. The integral 
equation then becomes a system of linear equations, which is solved by numerical 
methods. This method and several of its variants are related to Saltykov's method 
(1967) [55]. The way the data are grouped (the choice of the histograms bin-widths) is 
critical for such methods and may give unreliable results for moderate samples. 

One of the earliest analyses to be prompted by metallurgical problems is due to 
Scheil (1931) [56]. His method of tmfolding sphere size distributions is rather inflexible 
and involves lengthy calculations. ScheiTs solution of the integral equation that relates 
section distributions to sphere distribution, is based on a stepwise rather that on a 
smooth distribution function. The number of particles per unit volume is calculated for 
each size group, starting with the largest and proceeding down to the smallest size 
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group. The section diameters obtained from the plane of polish, with radii ranging from 
0 to Kmax (= Rmax), UTC divided into 15 size groups. Particles of one size group are 
considered to be of only one size. From the geometrical relationships and the general 
integral equation relating the section-size and particle-size distributions of spheres he 
derived the following expression 
Equation 2-53 

(A '. ),j = iK )j Ar 4(2]Y-(2i-if - + 

Considering the term in the square bracket in the above equation as Xjj Scheil’s basic 
equation can be written as 
Equation 2-54 

I -Mi 

ArX. . 

where each term to be subtracted is evaluated by 
Equation 2-55 

iff A.,., 

etc., in sequence from largest to smallest sizes. The weakest part in the Scheil’s 
treatment is in the estimation of (NyJj for j<15, the calculated values are required, and 
not the measured values {N^ ). . 

Schwartz (1934) improved on ScheiTs analysis in several respects [57]. The 
primary advantage of his method is that the number of particles in any size group can 
be calculated directly; i.e., there is no sequential dependence upon previous calculations 
of other particle sizes. Another improvement is that measured number of sections per 
unit area is used instead of the calculated values of number of particles per unit volume 
in the frnal equation. Schwartz didn’t use the general integral equation relating the 
section-size and particle-size distributions of spheres. Rather, the coefficients required 
for the calculation are deduced from the probability of a random plane intersecting a 
sphere at a certain section size. A disadvantage of Schwartz’s method is that the values 
of the coefficients change if a different number of intervals is selected. Hyam and 
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Nutting circumvented this problem by using a geometrical progression of size intervals 
[58]. 

Saltykov further developed Scheil’s method as modified by Schwartz. 
Saltykov’s method embodies the advantages of Schwartz’s modification, but at the 
same time is free of its shortcomings. A direct calculation of the number of particles in 
any size group may be made regardless of whether the number of particles in any other 
group is known or not. This method is discussed in detail in the next chapter. 

Saltykov’s area analysis (1967) attempts to improve on older methods and 
succeeds in several regards [55]. It does not require sequential calculations or tables of 
coefficients. It has the advantages of simplicity of calculation and finer differentiation 
of the small sizes. It is applicable not only to polydispersed systems of spherical 
particles, but in principle, to convex particles as well. In this analysis, all particles must 
have the same shape. In order to specify the size of a section, Saltykov adopted the ratio 
A! instead of the absolute area. Thus it is apparent that the correct value of is 

very important. A large number of particles must be examined on the plane of polish, 
especially if the particles have a complicated shape. The scale factor used to determine 
the class intervals for A! A^is based on a logarithmic scale of diameters with the 

factor 10'° *. The advantage of using A! A^to group the sections lies in the unique 


relationship of the particle shape to the //^versus yl//4^firequency distribution. By 

intersecting a given particle with many random planes, one can obtain a distribution of 
sectional areas which is characteristic of that particular particle. For spherical particles, 
this distribution can be calculated readily and the values are tabulated by Saltykov. For 
particles of other shapes, the calculation of the {N^) distribution versus A/ A^ would 

be extremely difficult. The computations follow along the lines described previously. 
The procedure is similar, except that one successively subtracts those sections 
contributed by larger spheres. The spherical particles have diameters Dj corresponding 

to [Ny )j . Saltykov derived the following general working formula which directly gives 

{Ny )j for any desired class interval. 


Equation 2-56 



ri.646l(jVj, -0.0415(;vJ,., 

L -o.onM.NX, -0 0079(ArJ„ -0.0038(Wj,_. -om^nX, 

' -0.0010(iVj,_, -0.0003(ArJ„ -0.0002(JV,),_,. - 0 . 0002 (iVj,„, 
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Note that j can have any integer value up to and including 12. For example, to 
calculate the value of {Ny , one would use just the first four terms in the brackets. It 

should be noted that for this particular derivation, the largest particle size corresponds 
to 7 = 1 the / = 7 , with subscripts varying through i - 1, / - 2 , i - 3, etc., down to 1 . 

The derivation by Spektor (1950) appears to be the first among the well-known 
analytical methods for determining particle-size distribution curves firom chord-length 
distributions [59]. His method requires the shorter computational time, but measuring 
and recording individual chord lengths is slower than applying a circular template for 
the area method. The one extremely sensitive point in both analyses is the measurement 
of the smallest particles. His final working equation apply to either a continuous or a 
discontinuous size distribution and is expressed as follows 
Equation 2-57 

) 4 r k); 

®\"[27-1 27 + 1 

where the {Ny )j and {n^ ). are the number of particles per unit volume and the number 
of chords per unit length of test line, respectively, for the indicated 7 -size. We can see 
form the above equation that no coefficients need be calculated for the (n^ )j terms. We 

merely select a convenient A value, then substitute it and the experimentally 
determined values of (n^ )j into the above equation. 

Lord and Willis (1951) also did the same type of analysis as that proposed by 
Spektor, in an attempt to obtain the volume fi:action and size distribution of spherical 
air bubbles in firost-resistant concrete [60]. But they pointed out that their method is 
applicable to any system containing “spherical dispersoids.” In treating the relationship 
between a distribution of chord lengths and a distribution of sphere diameters, they 
gave both an analytical solution and a graphical procedure. Their analytical solution is 
analogous to the Spektor’ s analysis, but their graphical solution is unique. The main 
steps in their graphical solution may be listed as follows: 

1) They show, for imiform spheres of diameter Dj , that the graph of the firequency 
distribution n^{l)o{ chord lengths / is a straight line from/ = 0 , ni(/) = 0 to a 
maximum at/ = Dj . Thus a monodispersed system of spheres is represented by 
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a triangle, the area of which equals the total number of chords per unit length of 
test line. 

2) For two monodispersed systems of spheres, there will be two triangles; for three 
systems, three triangles; etc. These triangles can be recombined, preserving the 
individual triangle areas, in such a way that the resulting figure represents the 
total system. 

3) A frequency distribution curve of chord lengths resulting fi-om penetration of a 
polydispersed system by a test line is constructed in the same maimer. The 
number of chords from each class interval contributes a triangular area, the 
dimensions of which are dependent on the diameter and number of spheres in 
the class interval. 

For a polydispersed system of spheres a complete chord length distribution is plotted as 
described above. Figure 2-1 illustrates such a curve, where the number of chords per 
unit traverse length divided by the interval length is plotted against the medians of the 
corresponding intervals. The intercepted distances along the vertical line AB are related 
to the number of spheres per unit volume in each length interval. They have shown by 
simple mathematical analysis that 
Equation 2-58 

^ 2(KN,), 

TV OA 

When the distance OA is set off 2/7t units on the / scale, the above equation reduces to 
Equation 2-59 

Taylor (1982) was the first to use a kernel estimator for estimating, by a smooth 
function, the density of the planar circular profiles which is then substituted in the 
integral equation to obtain an estimate of the density of the sphere radii [61]. Taylor’s 
method, which can be viewed as an operation smoother than the one produced by the 
histogram binning of Saltikov's method, relies also upon a good choice of the window 
width of the kernel smoother. Hall & Smith (1988) have shown that Taylor’s method 
can be efficient only if the number Na of observed planar circular sections is very large, 
which is rarely true in practical situations [62]. A kernel based approach for estimating 
the cumulative distribution of the squared three-dimensional radii has also been derived 
recently by Golubev & Levit (1998) [63]. Recently, Feuerverger & Hall (2000) studied 
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a thick-slice version of Wicksell's problem, and introduced an appropriate 
cross-validation criterion for selecting the kernel smoothing parameter [64], Nychka et 
al. (1984) used a cross-validated spline regularization method to estimate the density of 
the sphere radii [65]. The solution is obtained by approximating the unknown density 
by a smooth combination of natural splines, in such a way that its integral transform 
stays close to a histogram of the data, with a number of classes much larger than the 
one that should be employed to estimate optimally the density of the planar sections. 
Antoniadis (2001) developed a method based on wavelet transforms for unfolding 
sphere size distribution. The method consists of first converting the wicksell’s integral 
equation to a form suitable for application of thresholding wavelet methods and then 
deriving the asymptotic properties of their estimators [66]. 

The assumption of spherical grains is dubious for metal grains, which are 
known to be space filling. Scheil and Wurst (1936) determined the distribution of grain 
sizes based on actual grain shapes, rather than on grains having spherical or ellipsoidal 
shapes [67]. By successive sectioning of a sample of Armco iron, Scheil and Wurst 
obtained average section radii through 68 grains, firom which they calculated a 
statistical grain shape. They provided a table of coefficients based on 14 class intervals 
for this particular statistical grain shape. The derivation of the spatial size distribution 
firom the planar size distribution is similar to that described by Scheil for spherical 
particles. Paulaus (1963), who is interested in the magnetic and electrical properties of 
ferrites as influenced by their grain-size distributions, fixrther improved the Scheil and 
Wrast analysis. On the basis of comparison of average characteristics of real grains and 
cells, Paulus selected the pentagonal dodecahedron as the prototype shape for ferrite 
grains. Using Hull and Houk’s frequency distribution curve of diameters for the 
pentagonal dodecahedron, he calculated a table of coefficients in much the same way as 
Scheil and Wurst. Paulus offered both analytical and graphical methods [68]. The 
analytical method is based on Scheil’s original analysis, except that class intervals are 
spaced according to a geometrical progression. With the help of the latter 
simplification, plus the assumption of lognormal distribution, Paulus prepared graphical 
solutions for several special cases. The limitations inherent in a Scheil-type analysis are 
circumvented to some extent in Paulus’ graphical solution. 

Cahn and Pullman (1956) analyzed the lamellar structures, which is the 

extension of their procedure for spheres [69]. The method applies to structures such as 

pearlite, in which hard lamellae are embedded regularly in soft matrix. Instead of a 
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calculated or analytical distribution function of lamellae spacing, this method is based 
on a measured lamellae spacing distribution curve. Cahn and Pullman considered two 
experimental parameters: tme spacing cr, and random intercept spacing cr^ , the tme 

spacing is the perpendicular distance across two adjacent platelets, while the random 
spacing is the center-to center distance between hard lamellae along random test lines 
on the plane of polish. Using the geometrical probability principles and equivalence of 
volume ftaction and lineal fraction, they derived the following expression for volume 
fraction distribution function of the true spacing of lamellae depending on lineal 
distribution function: 


Equation 2-60 




diner. 


where g(cr ^ ) is the lineal fraction distribution function of random spacings. Similarly 


depending on number distribution function of random spacings they derived the 
following expression for volume fraction distribution function of true spacings: 
Equation 2-61 



= 3cr^nj.(a,)+cr^ 




In 1926 in Wicksell dealt with the case of ellipsoidal particles [70]. Wicksell chose 
the geometric means of the diameters of triaxial ellipsoids and of their profiles as the 
characteristic measure or “diameter”. He then found first that the frequency of profile 
“diameters” is related to the “diameter” of the central ellipse of the particle by the same 
fundamental Abel-type integral equation as he had found for spheres (1925). For 
prolate spheroids this was also the case if he chose the minor axis of both ellipsoids and 
section ellipses, or if for oblate spheroids he used the major axis. After rigorous 
mathematical analysis he concluded that exactly the same relations exist between the 
functions F (frequency of the central diameter of the ellipsoid) as in the case of 
spherical corpuscles. Hence applying the reductions deduced for the case of spherical 
corpuscles to the distribution of the minor diameters in section ellipses one can obtain 
the distribution of the minor diameters of the corpuscles, if they are all prolate 
ellipsoids of revolution. It is evident that the same result will hold if the corpuscles are 
all oblate ellipsoids of revolution, if we then apply the spherical reductions to the 


26 



distribution of the major diameters in the section ellipses. The result obtained will be 
the major diameters of the corpuscles. 

Till 1962, no further attempts seem to have been made to re-derive WickselTs 
solution for the size distribution of ellipsoidal particles. Dehoff (1962) modified the 
Schwartz-Saltykov analysis for spherical particles so as to obtain the size distribution of 
ellipsoidal particles of the same shape (axial ratio). This analysis is dealt in more detail 
in the next chapter. 

The difficulty with Wicksell’s and Dehoff s analyses of ellipsoidal particles is 
that it is not valid for a general population of ellipsoids because it was based on the 
rather strong assumption that size and shape were statistically independent. The general 
problem, which involves bivariate distributions of size and shape, has been solved by 
Cmz-Orive (1976) [71]. Spherical particles, or particles modeled by geometrical bodies 
of a constant shape (e.g. ellipsoids with constant axial ratios, cylinders with a constant 
height-to-diameter ratio, etc.), can be described by a single size parameter as proposed 
by Wicksell. The particulate phase is then described by a univariate size distribution, 
and the stereological problem consists simply in identifying this distribution from the 
univariate size distribution describing the set of profiles produced by a random plane 
section through the aggregate of particles. Particles exhibiting variation about a given 
type of shape as well as size variation, however, carmot be described by a single 
parameter. Such is the case of, ellipsoidal particles of variable size and eccentricities. 
Assuming that the spheroid centers are uniformly scattered within the specimen, the 
spheroids are isotropically and independently orientated about their centers and finally 
the size and shape of a randomly chosen spheroid are independent from its position 
within the specimen, Cruz-Orive solved the problem of estimating bivariate size-shape 
distribution of spheroids from the corresponding size-shape distribution of elliptical 
profile. This problem is explained in detail in the next section. 

A. M. Gokhale (1996) derived a fairly general and flexible statistical 

relationship between the apparent bivariate distribution of size and orientation of micro 

crack traces observed m metallographic sections (which can be experimentally 

measured), and the true bivariate size and orientation distribution of micro cracks in 

three-dimensional microstructure, without assximing anything concerning anisotropy of 

micro cracks, their sizes, orientations, thickness, and any correlations between size and 

orientation [72]. Prior to his derivation Seaman et al (1978) developed a numerical 

analysis algorithm to calculate the histogram representation of true bivariate size and 
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orientation distribution of micro cracks from measurements performed on a 
metallographic section [73]. But they assumed that the true distribution is rotationally 
symmetric with respect to a reference axis, i.e. there is an axis of symmetry. They 
formulated the problem in a discretized format and this discrete (rather than closed 
form analytical) treatment does not permit development of more efficient numerical 
algorithms. Gokahle’s formulation of the problem is not in discretized format. The 
apparent bivariate size and orientation distribution of micro cracks is related to the 
corresponding true bivariate distribution function through a double Abelian integral 
equation. This analytical treatment is convenient for calculation of the apparent 
bivariate size and orientation distribution of the micro crack traces for modeled tme 
bivariate distributions of micro cracks. Actually the micro cracks are modeled as penny 
shaped cylindrical objects with flat circular ends (platelet like). Hence this method is 
applicable to any system of circular disc shaped particles. The mathematical treatment 
of this method is discussed in detail in next section. 

The problem of unfolding bivariate distribution of unequiaxed particles 
modeled by triaxial ellipsoids of variable size and eccentricity parameters is 
indeterminate in its generality, because the elliptical profiles do not furnish sufficient 
information which permits a complete description of the ellipsoids [71]. Consider an a 
phase consisting of randomly oriented ellipsoids dilutely embedded in an opaque 
matrix. Let the principal semiaxes of an ellipsoid are denoted by a, b, c varying in the 
region (0<c:^^), (b:^<oo). Where, B is the largest value of b in the ensemble of 
ellipsoids. The two principal squared eccentricities are defined as X}^=l-(b/af and 
X2=l-(c/af' , so that 0< X]^< X2<1. Let the ellipsoids are described by the three 
dimensional random variable (c, xj^, x/j, consisting of a size component plus two 
shape components. Their joint probability density function g-(c, x/, X 2 ^) is defined in the 
region (0<c:^), (0<x/<X2^<l). The population of random, plane elliptical profiles is 
described by the bidimensional variable (m,)/^), with joint probability density function 
f(m,y^) defined in the region (0<m^), (0<y^<l). Where, m is the minor axis of the 
profile andy its squared eccentricity parameter. The general version of Equation 3-29 
can be expressed as. 
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Equation 2-62 



Where f denotes the conditional joint probability density function of (m,y^) on the 
ellipses produced by a random plane T in a fixed ellipsoid K of parameter (c, x/, X 2 ^), 
whereas ^*=gfc, x/, X 2 ^\TfK). 

Let us suppose that all the ellipsoids in the a-phase are spheroids. Hence, the 
whole probability mass is concentrated in the two plane rectangular regions {(0<c:^), 
(0<xi^= X2^<1)} and {(0<c:^), (0= xi^<X 2 ^<l)}. In the former case, g takes the form 
gp(c, x^, x^)^p(c, x^), which describes the subpopulation of prolate spheroids. In the 
latter case g becomes, go(c, 0, x?) =go(c, x^), which describes the oblate subpopulation. 
Clearly, gp and go must satisfy the identity 
Equation 2-63 

0 0 

For the mixture of spheroids as the one described above, the general expression 
Equation 2-62 reduces to 
Equation 2-64 

/(m, y ' ) = J /<, j 

= j foSl^ + J fWpdc 

M in 

Clearly, when go=0 (all spheroids are prolate), the equation reduces to the equation for 
prolate ellipsoids, and when gp=0 (all the spheroids are oblate), the equation reduces to 
the equation for oblate spheroids. It is shown in Equation 3-27 and Equation 3-29, that 
these equations have unique solution. However, Equation 2-64does not admit unique 
solution for and gp. To check this let us p\itf=qf+(l-q)f, where ^ is a fixed but 
arbitrary number between 0 and 1 and consider the following equations. 
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Equation 2-65 


of = \foS’oda 

M 

il-g)f= jdx^ jf^gldc 
/ « 

The above equations are of the form Equation 3-49 and go and gp can be solved for 
unique solution as per Equation 3-5 1 . It is evident that the two functions so obtained 
satisfy Equation 2-64, and hence they are solutions of it. However, one can obtain as 
many pairs (go, gp) satisfying Equation 2-64 as one pleases, by simply choosing 
different values of q. It is clear from the above discussion that the problem of “mixture 
of spheroids” and hence the problem of “triaxial ellipsoids” is indeterminate. 

Warren and Durand proposed a computer-modeled approach, which they named 
as “guess-and-compare process” to determine the 3-D particle size distribution in 
systems of uniform, identifiable particle shape. The methodology of this process is, 
guessing a 3-D size distribution and then, by computer modeling producing a 2-D 
result that can be compared with the experimental distribution. This process of guess 
and compare is repeated until satisfactory agreement is achieved. 

In case of polyhedral shaped particles, the shqje of the profiles will be triangle 
or any other polygon depending upon the order of the polyhedra of the particles. In case 
of cubes profiles will be anything between triangles and hexagons. So it does not make 
any sense to define any particular “diameter” in this case. The problem of estimation of 
the size distribution of cubical particles from the distribution of section profiles 
remained undetermined till J.Osher and M. Nippe (1997) proposed a method exploiting 
size and shape of the observed profiles [74]. They actually derived a methodology that 
is applicable to any homothetic particle system. As mentioned earlier the section 
profiles no longer form a system of homothetic figures, hence they carry additional 
information than can be expressed by one size parameter. They made two observations 
on planar sections: 1) measuring the maximum breadth of the profiles (caliper 
diameter) and 2) counting the number of vertices of the profiles. They also have shown 
that the size of a single cubic particle can be determined immediately from the shape of 
a given section profile if the number of vertices of that profile is greater than 3. 

All the stereological unfolding procedures discussed above are model-based 
methods, i.e. the objects are assumed to have a particular size, shape, and orientation or 
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the objects are modeled (e.g., sphere, ellipsoid, cube, etc.). The data for all the above- 
mentioned unfolding methods comes fix)m the 2-D measurements, which always 
contain sampling errors and bias. In almost all methods of deconvolution, calculation of 
the 3-D particle size distribution from 2-D profile distribution essentially involves a 
numerical solution to an Abelian integral equation. These numerical calculations still 
amplify the already present sampling errors and bias in the 2-D data and hence the final 
solution consists of more errors than in the raw data. 

The other category of stereological methods, design based methods, are 
characterized by counting mles and sampling schemes that have been designed to 
eliminate considerations of the size, shape, and orientation of cells being counted. 
Design based techniques for counting requires not a priori knowledge about the size, 
shape or orientation of the objects being counted. Serial sectioning is one of the well 
known design based techniques. Conventional serial sectioning has been successfully 
used for quantitative characterization of microstructures, viz. to study the 3-D 
morphologies of the features that appear on microstructure, to estimate the number of 
particles per unit volume of space and to quantify the topological connectivity of 
microstructural features in 3-D. Conventional serial sectioning involves digitally 
recording a microstructural field of view, removing a small thickness of the specimen 
by controlled polishing, again recording the field of view exactly below the first one on 
the new section. This procedure is repeated to generate a stack of large number of 
aligned serial sections. The 3-D information of the microstructural features can be then 
obtained from this stack of serial sections. Since this technique involves serious edge 
effects and sampling bias, it has not been used for the estimation of 3-D particle size 
distribution. In order to get accurate measurements of particles sizes, microstructure has 
to be observed at sufficiently high magnification and the distance between consecutive 
serial sections has to be significantly less than the average size of the particles/grains 
under observation. At high magnification the field of view hardly contains 25-50 
particles/grains. In practice at the most about 50-100 serial sections are generated. 
Therefore, the microstructural volume generated by conventional serial sectioning 
technique contains, at the most, a hundred particles or so. Such small statistical sample 
size is not sufficient for reliable estimation of the complete particle/grain size 
distribution. The edge effect also brings a lot of systematic bias into the raw data. This 
is because of the fact that larger particles are more likely to intersect the jframe 
boundaries. 
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. A. Tewari and A. M. Gokhale (2000) developed a method by the application of 
montage-based efficient serial sectioning technique for estimation of particle size 
distribution, that overcome the shortcomings of conventional serial sectioning 
technique [75]. The basic idea behind this newly developed technique is that, first a 
small microstructural volume is reconstructed from a stack of conventional serial 
sections and then many contiguous small volumes are reconstructed surrounding the 
initial volume and finally all of them are pasted together after perfectly matching their 
boundaries to generate a large (about 25-100 times larger than the volume developed by 
conventional technique) three dimensional microstructural volume at high resolution 
from the same number (50-100) of serial sections as that in the conventional serial 
sectioning. The beauty of this new technique is with the same metallographic effort as 
in the case of conventional technique, one can generate a microstructural volume 
containing few thousands of particles. Such a large statistical sample is sufficient for 
reliable, unbiased, assumption-free and efficient estimation of 3-D particle size 
distribution. 
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Figure 2-1 Distribution curve of chord lengths constructed by graphical method [18] 
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CHAPTER 3 


COMPUTATIONAL BACKGROUND 

3.1 SIZE DISTRIBUTION OF SPHERES 

3.1.1 INTRODUCTION 

Before going into detailed discussion about unfolding the size distribution of 
spheres from measured size distribution of circular profiles, a brief note is given here 
on how the profile sizes are distributed when a single sphere of radius R is randomly 
sectioned many times [76]. 

Consider a sphere of radius R, which is embedded in a cube of matrix material 
and is placed in x, y, z coordinate system with its center at the origin (See Figure 3-1). 
Now, a large number of equidistant section planes, which are parallel to x-y plane are 
placed into this cube of material. Only those sections, which lie between the upper and 
lower parallel tangent planes will cut the sphere inside the cube. It is evident that the 
number of section planes that cut the sphere is proportional to the diameter of the 
sphere, i.e. 2R (in this case). Each of these section planes will produce a circular 
profile. The radii of the circular profiles thus produced are a fimction of the position of 
the section planes with respect to the z-axis, as shown below. 

Equation 3-1 

r(z) = ^(R^-z^) 

The size distribution of profiles is characterized by determining the probability 
density with which a profile radius p falls into the range 
Equation 3-2 

r>ri >{r-dr) 

Since the sectioning planes are parallel to x-y plane, the probability measure must be 
related to the distance along the z-axis between the plane that will produce a profile of 
radius r and that, which produces a profile of radius (r-dr). Let us denote this distance 
by dz. From the rules of geometric probability, the probability of sectioning any object 
with a plane is the ratio of the tangent diameter of the object to the tangent diameter of 
the containing space. Hence, the probability of having q in the class {r, (r-dr)} is given 

by. 
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Equation 3-3 


Pr{r > r, > (r — dr)\ — 


dz! R for r < i?l 
0 for r>i?J 


It should be noted that the slice thickness dz is not constant, but is a function of r as 
shown in Figure 3-2. From the definition of dz one can write. 

Equation 3-4 

dz = z{r - dr) — z(r) 


From the equation of the circle we can derive, 

Equation 3-5 

z(r) = -^{r^ 

Equation 3-6 

zir-dr) = ,l{R^ -{r-drf) 

Inserting Equation 3-5 and Equation 3-6 in Equation 3-4, we get. 
Equation 3-7 

dz = ^[R^-{r- drf ) - - r' ) 

Differentiating Equation 3-7 with respect to r, we get. 

Equation 3-8 

dz _ —r 
dr 


Inserting Equation 3-7 in Equation 3-3 the probability density becomes. 
Equation 3-9 


Pr{r >r^ > (r - dr) | i?} = — -(r-dr)^] - ~r^) 

R L 


Inserting Equation 3-8 in Equation 3-3 the probability density becomes. 
Equation 3-10 


^)(r I R)dr = 



dr 


In the above equation the sign of the right hand side term is changed to give positive 
probabilities. In Equation 3-9 and Equation 3-10 the argument of the left hand side 
function is extended to note that it depends on both the sphere size and the profile size. 
It should be noted that Equation 3-10 is continuous version of Equation 3-9 and it 
becomes Equation 3-9 if ^(r | R) is integrated between the class limits {r, r-dr}. 
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Equation 3-11 


Pr{r, r - c?r} = J^(r | /?)c/r = Pr{r | i?} 

r-dr 

The form of the curve <j>{r\R) is shown in fig (6.3a). Since (f){r) is a probability 
density fimction, if Equation 3-10 is integrated over the range 0 to R, the result must be 
unity. A more realistic representation of this function is perhaps a histogram as shown 
is fig (6.3b) whereby dr may be chosen as fine as desired ( Jr -> 0 at the limit). The 
values Pr{r\R} of the histogram columns is obtained by Equation 3-9. 

From the above analysis we find that, although all profiles are equally likely to 
occur, the probability density of profile radii varies with profile size, i.e. larger profiles 
are more likely to be found on section than smaller ones. This is not paradoxical. 
Because, in defining the probability of obtaining a certain profile we proceed at 
sectioning the sphere at equal intervals along the z-axis normal to the sectioning plane. 
The increments dr(z) in profile radius will then become smaller as the section plane 
moves down as shown in Figure 3-3. In defining the probability density of measuring a 
certain profile radius we assume equal increments of r by dr. This produces larger 
intervals dz(r) between the slice boimdaries as r becomes larger as shown in Figure 3-4. 
In the above analysis we have assumed that one sphere was sectioned many times. This 
assumption is similar to assuming a monodispersed system of spheres being randomly 
embedded in a large block and this block being sectioned once, or a few times at great 
distances. 

3.1.2 METHODOLOGY 

3.1.2.1 DISCRETE MODEL 

Consider a polydispersed system of spheres, i.e. a system of spherical particles 
of varying size, such as the one shown in Figure 3-5. If this population of spheres is 
randomly sectioned with a plane, one gets a distribution of circular profiles. There are 
basically three factors that may influence in some way or another the number of 
profiles contributed to each profile class. They are, 

(a) The numerical density of spheres occurring in each class of the sphere size 
distribution; 

(b) The probability of hitting a sphere of a certain size with a random sectioning 
plane, which is related to its (mean tangent) diameter; 
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(c) The frequency distribution of profiles derived from each clziss of spheres. 

Let us assume that the sphere and profile size classes have the same interval, i.e. 
AR = Ar. Hence, the sphere and profile radii can be expressed as integer multiples of 
/!/? or zlr, as follows. 

Equation 3-12 

Rj=j.AR = j.Ar 
Ki = LAR = /.Ar 


Let the number of spheres of size Rj per unit volume of the system be (Ny)j and the 
number of section profiles that fall in the size class / be (7/^,-. Let the radius of the 
largest sphere is Rk. From the basic stereological principle regarding particle size 
distributions Equation 2-48 we can write. 

Equation 3-13 


i 



Where the subscript j refers to the sphere size and the subscript / refers to the section 
size. The term ^ {N^ ),. y expresses the number of sections of all sizes per unit area 


obtained from particles of a fixed size J. The contribution of spheres of class j to 
profiles of size class / can be determined by the profile size frequency given in 
Equation 3-9, as follows, 

Equation 3-14 

i 



It should be noted that not all sphere classes may contribute to {Na):. It is evident that 
spheres whose radius Rj < cannot give rise to any profiles of size class i. Keeping 

this in mind it can be written as. 

Equation 3-15 

(NA-t(NX 

Jm 

Inserting Equation 3-14 for in Equation 3-15 we get. 
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Equation 3-16 


{NA='Li‘^.)r2 

jmi 

Dividing the above equation with Na on both sides. 
Equation 3-17 



Where f(ri) is the frequency of observed profiles in size class i. 

Now considering the basic stereological relation according to Equation 2-48 
Equation 3-18 

N^=Ny.2R 

Where R is the mean radius of all spheres, which is given by 
Equation 3-19 

Jy/y 

Inserting Equation 3-18 for Na in Equation 3-17 we get 
Equation 3-20 




Where F(Rj) is the frequency of spheres of size class Rj. Inserting Equation 3-12 for Rj 
and n in Equation 3-20 we arrive at the following equation relating the frequency 
distribution of profile sizes in 2D and the frequency distribution of corresponding 
sphere sizes in 3D. 

Equation 3-21 

/(o=ii;p 

J M 
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Now, the problem is to work back from the above equation, i.e. to estimate the 
frequency distribution of sphere sizes in 3D from the measured frequency distribution 
of profile sizes on 2D section plane. 

It should be noted from Figure 3-6, that the circular sections have a distribution 
of sizes, even though they are all obtained from the spheres of uniform diameter. This 
is due to the fact that the sectioning plane intersects each of these spheres at different 
distances from their centers, giving rise to sections of different sizes. Another important 
fact is that in a polydispersed system of spheres, the spheres of different size will 
contribute sections of same size. In order to find the distribution of three-dimensional 
spheres from the distribution of two-dimensional sections, the source of each of the 
sections must be determined. From Equation 3-14 we can write the following general 
expression 
Equation 3-22 

from which {Ny )^. may be determined. 

Let us assume that the polydispersed system of spheres under consideration is 
separated into k class intervals. The number of spheres per unit volume in each class 
interval is 

(Ny\,(Ny\, .(Ny\ 

and their corresponding radii are 

ri,r^ ,r. 

The section profiles observed on the plane of polish are also separated into the same 
number of (k) of class intervals. The number of circular sections per unit area in each 
class interval is 

-My 

where the section radii vary from 

Oro r,,r, tor^, ,r*_i to 

respectively. Figure 3-7 illustrates the possibilities arising from a random plane 
intersecting a polydispersed system of spheres with kr=5. It can be clearly seen fix>m 
Figure 3-7 that spheres of size r, yield sections of radii Keeping this fact in 

mind we can write 
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Equation 3-23 


Jm 

In order to use Equation 3-22 to calculate the number of spheres per unit volume in a 
specific size class, we have to determine the number of sections having the same radius 
as their spheres, i.e. • Since the sections in the last class 

interval are contributed only by the spheres of largest size we start with the largest 
particles and their sections and proceed to the smallest size. . From the above 
discussion and from Figure 3-7 we can easily write the following equations. 

Equation 3-24 

ANA.,-(NA-tj. -(NA-^J.-, 

(N , = [N, ), - (N, )„ - (N, X.,, - - (AT J,3 - (JV, X, 

In the above set of equations, the first term on the right hand side of each equation is a 
measured quantity. Since {N^)^ is a measured quantity, is known. We can 

calculate Pr{r^ 1 } from Equation 3-9. Substituting this probability value, {N^ and 

the measured in Equation 3-22, we can calculate {Ny)^ , i.e. the number of three- 
dimensional spheres per unit volume of the largest class interval. Having the value of 
in hand one can now proceed to calculate {Ny )j^_, . The probability of hitting a 

size-k sphere to get a size- (k-1) section can be calculated using Equation 3-9. 
Using this value in Equation 3-22, the value of (A^^ )*_, * can be calculated, which when 
subtracted from the measured value yields (as shown in Equation 

3-24). The probability for hitting a size- (k-1) sphere to get a size- (k-1) section 
( Pk-i k-x ) can be calculated using the Equation 3-9. Using this probability and knowing 

A )*-i k-i ’ )t-i be easily calculated from Equation 3-22. 
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By an extension of this procedure, step by step, until all class intervals are 
exhausted, we can calculate the number of spheres per unit volume in all the size 
classes. Now the frequencies of sphere size distribution can be easily calculated from. 
Equation 3-25 

L 

/»! 

It can be clearly seen that in the calculation of {Ny )j each subtracted term depends on 

the previously calculated values of etc. This leads to the 

accumulation of errors towards the smallest particle size and the calculated number of 
smallest particles may be inaccurate, especially if the number of size classes is more. 



3. 1.2.2 CONTINUOUS MODEL 


In the above model the population of spheres is considered to be made of 
discrete size classes and the profiles are sorted into likewise discrete profile size 
classes. For analytical purposes it is convenient to allow both the sphere size 
distribution F(R) and the profile size distribution f(r) to be continuous by letting both 
dR and dr become infinitesimally small, i.e. dR-^0,dr-^0 at the limit. From 
Equation 3-10 and Equation 3-1 1 we can write. 

Equation 3-26 




r.dr 


If we now replace the summation in Equation 3-21 with an integral and considering 
Equation 3-26 we obtain the following continuous distribution. 

Equation 3-27 

f(r)dr = i f“ F(r). ,, -rJR 

Dividing the above equation on both sides with dr we obtain the following fimdamental 
probability density fimction for profile size, relating the distribution of spheres, and that 

of profiles. 

Equation 3-28 
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Allowing for continuous distribution of sphere sizes, the above equation can be 
rewritten in its final form as 


f(r)= ^ F{rM’-\R)‘‘R 

It should be noted that the upper limit of integration is changed to oc instead of Rmnx to 
allow for continuous distribution of sphere size. 



3 . 2 SIZE SHAPE DISTRIBUTION OF VARIABLE SPHEROIDS 
3.2.1 INTRODUCTION 

The constraints of the surroundings often impose a particular shape on cells and 
their constituents, especially during their growth stage. Thus cells, which are basically 
spherical in shape, become polyhedral when they are densely packed, and the nuclei 
simulate the shape of the cell, becoming flattened or oblate ellipsoids in flattened cells, 
or prolate ellipsoids in columnar or elongated cells. There is a certain error introduced 
if one substitutes, a “volume-equivalent” sphere for non-spherical particles, but this 
error is predictable. In case of higher order polyhedra, as one can expect, this error is 
rather small, but for shapes such as cubes, octahedra or cylinders it is appreciable. With 
ellipsoids the error obviously increases as their eccentricity (or axial ratio) deviates 
from 1. 

Particles of spherical shape or those modeled by geometrical bodies of a 
constant shape (e.g. ellipsoids with constant axial ratios, cylinders with a constant 
height-to-diameter ratio, etc.), can be described by a single size parameter. The 
particulate phase is then described by a univariate size distribution. In this case the 
stereological unfolding problem consists simply in estimating the univariate size 
distribution of particles from the univariate size distribution of profiles produced by a 
random plane section through the aggregate of particles. Very often, this assumption is 
unrealistic, and complete information relative to the joint variation of size and shape is 
called for. In many cases particles exhibit a variation about a given type of shape as 
well as size. Such a system of particles can’t be described by a single parameter. Such 
is the case of ellipsoidal particles of variable size and eccentricities. If one allows the 
size and the shape of the ellipsoidal particles to vary, one will have to establish a 
relationship between a particle size and shape distribution and the corresponding profile 
size and shape distribution resulting from random sectioning. This matter is not a 
simple straightforward one. This is because of the fact that even if all the particles are 
of same shape, the profiles resulting from random sectioning will not have the same 
shape. For example, consider a prolate ellipsoid of axial ratio 2:1. When this particle is 
randomly cut with a plane, the resulting profiles will have axial ratios ranging from 1:1 
(circles) to 2:1 (ellipses), depending on the spatial orientation between the particle and 
the section plane as shown in Figure 3-8. 
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The problem of unfolding a bivariate distribution describing a phase of variable 
ellipsoids of revolution (spheroids) from the corresponding bivariate profile distribution 
is determinate provided the spheroids are all of the same type, namely, either all prolate 
or all oblate (See Figure 3-9). This is because of the fact that the random, elliptical 
profiles can fiunish complete information only if the ellipsoids are spheroids of the 
same family. Prolate ellipsoids are generated by ellipses revolving around their major 
principal axis, whereas oblate ellipsoids are generated by ellipses revolving around 
their minor principal axis. If the particulate phase consists of a mixture of prolate and 
oblate ellipsoids, then the relative proportion of each type of spheroids in the mixture 
introduces an extra parameter, making the problem indeterminate. Modeling a phase of 
particles of variable size and shape by variable spheroids of the same type cut by a 
sectioning plane seems appropriate for three main reasons. (1) Often, the model will be 
much more realistic than a sphere’s model, or a model consisting of spheroids of 
identical shape. (2) The profiles are always of the same type, namely ellipses, thus 
making the problem manageable in theory and practice. (3) The stereological problem 
is determinate (If one uses line probes instead of plane probes, the problem becomes 
indeterminate, because intercept length distributions are essentially univariate.) A few 
assumptions are made in setting up the model. (1) The spheroid centers are uniformly 
scattered in the specimen. (2) The spheroids are isotropically and independently 
orientated about their centers (3) The size and shape of a randomly chosen spheroid are 
independent from its position within the specimen. All the foregoing assumptions will 
be fulfilled only approximately, because the spheroids cannot overlap. 

3.2.2 METHODOLOGY 

3.2.2. 1 CONTINUOUS MODEL 


Any spheroid (either oblate or prolate) is fully determined by giving its major 
and minor principal semi axes a and i respectively. Here in this model both a and Z? 
are variable [71]. For both types of spheroids the eccentricity parameter is defined as 

1 _ _ ^ so that 0 < < 1 . To make the manipulations simpler a prolate spheroid is 

U; ’ 

described via the bidimensional random variable {b,x^) and an oblate spheroid via 
(a,x ^ ). Clearly, a or 6 is a size variable, whereas x^ (which is the square eccentricity 
parameter) is a genuine shape variable. The respective joint probability density 
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functions g[b,x^) and g{a,x^), exist and are defined in the rectangles 
^ < a < B,0 < x^ < l| and {o < 6 < 5,0 < < l| respectively. Random plane probes 

(random with respect to the particles, but not necessarily with respect to the specimen) 
determine elliptical profiles in the transected spheroids. The elliptical profiles are 
described by the variable with joint probability density function 

defined in a rectangle |o < m < 5,0 < y^ < l| if the parent spheroid is prolate, and by 
the variable with joint probability density function /(M,y^) defined in a 
rectangle {o < M < 5,0 < < l| if the parent spheroid is oblate. Where M and m 

represent the major and minor principal semi axes of the elliptical profile and 


.y‘ 


: 1 - 


\Mj 


is its square eccentricity parameter. 


From the conditional probability calculus, for prolate spheroids one can write the 
following expression. 

Equation 3-29 


B 1 

/(m,y^)= \db\f[my \b,x\T K)g{b,x^ \T ^ K)ix^ 


In the above integrand the first function represents the conditional probability density 
function of ) on the set of ellipses produced ( i.e. probability of getting an ellipse 
of size and shape(/w,>’^) ) when a prolate spheroid K of fixed size b and shape x^ is 
transected by a plane T. The event of plane T transecting the spheroid K is denoted by 
T'\ K . The second function in the above integrand represents the conditional 
probability density function of (6,x^) on the population of transected spheroids (i.e. 
probability of T hitting the spheroid K of size and shape (6,x^)). In case of oblate 
spheroids Equation 3-29 becomes, . 

Equation 3-30 

/(W.y’)= ]db I t K)g{a,x‘ \ T t k)1x^ 

M 

The functions in the above integrand have the similar meaning as explained above in 
case of prolate spheroids. 

Consider a fixed spheroid K, with its center at the origin of coordinates, and its 
rotatory axis along the Z-axis. When a plane T , defined by the spherical polar 
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coordinates {p,<p,d) (of the foot of the perpendicular drawn from origin to T) 
intersects K, the size and shape of the resulting elliptical profile have the following 
expressions; 

If AT is a prolate spheroid. 

Equation 3-31 


m = 




=x^sin^ 9 


If AT is an oblate spheroid, 

Equation 3-32 

2 _ sin^ 9 
1-r^cos^^ 

Since the element of invariant measure for random planes in three dimensions 
equals s\xi9 dp d(f> d9 [1], the joint probability density fimction of (p,^) conditional on 
the event that the random plane T hits AT, is proportional to sin9 . For prolate AT it 
becomes. 

Equation 3-33 

’ t[p,9\b,x\T^ k)= [lEg [/f (0)F‘ sin 9 , 

For oblate AT, 

Equation 3-34 


M = 


a 


\~x^ COS" 9 


for p e [- H{9\H(d)\9 e [ 0 , 7 t I2\andt{p,9 i b,x\T t A:) = 0 otherwise. ^"Equation 
3-33 and Equation 3-34 the fimction H(9), represents half the caliper diameter of the 
spheroid in the direction 9, and Eg{H{9"^ is its mean value over all the directions. For 

prolate spheroid half the caliper diameter is given by, 

Equation 3-35 

H(9) = b{l-x^Y^{l-x^s\n^9Y^ 

For oblate Sf^heroid it is given by, ■ 

Equation 3-36 

Af((9) = a(l-x^cos^^)^ 


46 



The mean caliper diameter of a spheroid is given by, 
Equation 3-37 


’‘A 

^e[H{0)]= jH{0)sinede 

0 

After making the change of variable Equation 3-31 in Equation 3-33, the first 
probability density function in the integrand of Equation 3-29 becomes: 

Equation 3-38 

f(m,y^\b,x\TrK)=[2E,{H)]- 

for m € € [o,;c^ j and / = 0 otherwise. For oblate spheroids making the change 

of variable Equation 3-32 in Equation 3-33 the first probability density function in the 
integrand of Equation 3-30 becomes: 

Equation 3-39 


m 

1 

1-y' 


1 

X 

1 

1 

r4 

1 



a,x\T^ k)=[2E,(h)Y 


M 


^ 1 


1 

ts 

1 

1 

I 


y 


for M e (0,a],>'^ e [o,jc^ j and f = 0 otherwise. One can notice from Equation 3-38 
and Equation 3-39 that these two functions factorize as product of two independent 
functions as f^{m,orM )- It can also be noticed that the bounds of m, (or M), 

do not involve y^, and vice versa. This indicates that for a fixed spheroid, the size {m, or 
M) and shape (y^) of an elliptical profile are independently distributed, whatever the 
type, size and shape of the parent spheroid. 

Considering the probability that the plane probe T hits the spheroid K is 
proportional to half the mean caliper diameter of the spheroid, the second probability 
density function in the integrand of Equation 3-29 can be expressed as: 

Equation 3-40 

g{by I r t Ji:)= 


for 6€(0,5 ],x^ €[0,l) and g-^O, otherwise. For oblate spheroid the second 
probability density function in the integrand of Equation 3-30 becomes: 

Equation 3-41 

g{a.x^\T'tK)^ [e,E, (ff )]-' [E, {H)]g(a,x ' ), 
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for a € (0,5],-ii: € [0,l) and g = 0 , otherwise. In Equation 3-40 and Equation 3-41 the 

function E^Eq{H^—{h^ denotes half the mean of mean caliper lengths in the 

population of spheroids. It is used as normalizing constant in the above two equations. 

Inserting Equation 3-38 and Equation 3-40 in Equation 3-29 and Equation 3-39 
and Equation 3-41 in Equation 3-30 the integral equations expressing the profile 
densities in terms of particle densities become as follows: 

For pro late spheroids. 

Equation 3-42 



for m e (O,^],;;^ € [0,l),and f = Q otherwise. 
For oblate spheroids. 

Equation 3-43 



for M € € [0,l),and / = 0 otherwise. It can be noted fi-om both Equation 3-41 

and Equation 3-42 that the first integral depends exclusively on particle and profile 
size, whereas the second integral contains both the bivariate size-shape distribution and 
the shape parameters. 

In order to unfold the three-dimensional particle size and shape distributions 
firom two-dimensional profile size and shape distributions one needs the inverse 
relations of Equation 3-41 and Equation 3-42. Both of these equations are double 
integral equations of Abel-type. In order to transform these two equations into a single 
one of a standard type, the following substitutions are made: 

For prolate spheroids. 

Equation 3-44 

and g{p,x^)=r(b,x) 

For oblate spheroids. 

Equation 3-45 

AS 


After making the above substitutions in Equation 3-41 and Equation 3-42, the final 
integral equation whose solution is sought becomes, 

Equation 3-46 


H J /7i 7J ’ 


r{s,t) 


which can be regarded as a double integral equation of Abel type. To solve this 
equation put the second integral as. 

Equation 3-47 

V V/' -V* 

Equation 3-47 is an ordinary Abelian integral equation. Sneddon in 1966 [77] has 
solved this equation and the solution is 
Equation 3-48 


^ ^ rrat J /..2 .2 


TV dt j 

On the other hand. Equation 3-46 after inserting Equation 3-47 in it becomes, 

Equation 3-49 

Equation 3-49 is in the same form as Equation 3-47, hence its solution is 
Equation 3-50 

Inserting Equation 3-50 for the function i{/(s,v), in Equation 3-48, we obtain the 
following unique solution of Equation 3-46 
Equation 3-51 

AH d 




f'f 

' d ^ 





By using the following mathematical relation 
jF(w,5)is = 


ds 
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repeatedly, the right hand side of Equation 3-51 can be expanded in a form suitable for 
our purposes as. 

Equation 3-52 



If the function (f) satisfies the following conditions 

(p{B, v) s 0 and ^(u,l) s 0 for all u,v 

Equation 3-52 reduces to 
Equation 3-53 



Taking Equation 3-44 into account the above mentioned conditions are equivalent to 
Equation 3-54 

/(5,y^)s0 and lim =0, 

respectively. The equalities in Equation 3-54 imply that the population hardly contains 
big (size == B), very unequiaxed (y^ ->1) profiles. This assumption is a realistic one. 
Under the above mentioned conditions firom Equation 3-44 and Equation 3-53, the 
solution of Equation 3-42, for the system of prolate spheroids becomes 
Equation 3-55 



for & e (0,5],x^ g [0,l), and g s 0 otherwise. 

In case of oblate spheroids, if the joint size-shape probability density function of 
elliptical profiles satisfies the following expression. 

Equation 3-56 

/(5, y ^ ) * 0, and lim (l -y^)^ f{M,y^ )J= 0 
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then fix)ni Equation 3-45 and Equation 3-53 the solution of Equation 3-43 is given by, 
Equation 3-57 



for a e e [0,l), and s 0 otherwise. 

3.2.2.2 DISCRETE MODEL 

In practice, the profile size and shape distributions are obtained by distributing 
the measurements into size and shape classes. The resulting distributions are discrete 
(See Figure 3-10) in nature. A similar procedure as that described in case of spherical 
particles can be derived for unfolding discrete size-shape distribution of spheroids. 

The size variable b ox a , of prolate or oblate spheroid is assumed to vary in the 
range (0, B), where B is a constant bigger than, or equal to the largest value of b (or of 
a). This range is divided into s classes of equal width A = B/s, (5>l). The shape 

component varies in the range (0 ,1), which is divided into k classes of equal width 
1/k. Thus the rectangular domain of variation of {b,x^) or {a,x^) is divided into a grid 
comprising of sxk classes, each class being represented by a rectangle of sides A and 
1/k (See Figure 3-1 Ob). It should be noted that each of the spheroids in the particulate 
system under consideration belongs to exactly one of such classes. Any spheroid that 

satisfies the inequalities (i - 1)A < b{or a) < zA and ^ ^ , where i = 1, 2, 

..., s and j = 1, 2, ... , A:, belongs to ij th class and is called an zy-spheroid. The number 
of zy-spheroids per unit voltime of specimen is denoted by Ny{i,j), so that 

j * 

^ ^ y ) = JVy , which is the overall numerical density of spheroids. The elliptical 

M y-1 

profiles are also classified by means of the same size-shape grid used for the spheroids 
(See Figure 3- 10a). Thus any elliptical profile that satisfies the inequalities 

(a -1)A < m(orM) ^ oA and < y' < ^, where a= 1, 2, ..., s, and 1 ,2 , 

... , k, belongs to aj0 th class and is called a^ellipse. The number of a/B^Uipses 
observed on random section per umt area of section is denoted, by (a, that 

51 mm fin 


141936 ' 


Xl-S ^ which is the total number of elliptical profiles per unit area of 

a«“l 

section. From basic stereological relationships and from geometrical probability theory 
the ellipse number densities are related to the spheroid number densities as follows: 
Equation 3-58 

i«a 


for each a~ 1, 2, ... , /?= 1, 2, Clearly this relationship is the discrete form of 

the continuous equations Equation 3-55 and Equation 3-57. It can be seen that this 
equation is analogous to Equation 3-22 for spheres. Here in the above equation the 
coefficients are size-correctors and they are independent from A- as in the case of 

spheres- and also from the spheroid type. They can be calculated from the following 
expressions [78]: 


Equation 3-59 


7/ . ^ 


V 1 ^ 


1 

I — 

V 2y 




(. 3 ^ 

I — 

4j 


, for a <i 

for a = i 
for a >i 


It can be noted that this equation is similar to the one given for spheres Equation 
3-9. In Equation 3-58 the coefficients qjp are shape correctors, and as one can expect 

they obviously depend on number of shape classes ‘A:’ and also on spheroid type. These 
coefficients can be calculated from the following expressions: 

For prolate spheroids. 
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Equation 3-3 




for p<j 



for P = j 


0, 

V 

for p>j 

where fif) = arc tanh 

-2y0 + l)^ 

For oblate spheroids. 

Equation 3-4 


for P<j 



for p = j 


0, 

for p> j 


In order to estimate the size shape distribution of spheroids from the measured 
size shape distributions of elliptical profiles one has to reverse Equation 3-1. The 
reverse relation reads as follows; 

Equation 3-5 

j) = A-' Y. A («. 

a=i J5=j 

for all i = 1 , 2, . . . , A-; 7 = 1 , 2, . . . , yt. In the above equation the coefficient is the m th 
element of the inverse of a matrix whose cd th element is just ; likewise, the 
coefficient is the pj th element of a matrix having as the jp th element. It 

should be noted that the density of spheroids in any class can be directly obtained using 
Equation 3-5, without any need to calculate the density of spheroids in the preceding 
classes. The main drawback of this approach is that it might require very large profile 
samples. 
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3.3 SIZE AND ORIENTATION DISTRIBUTION OF CIRCULAR 
DISKS 

3.3.1 INTRODUCTION 

In all the above analyses it was assumed that the particles have random 
orientations (in other words the structures are isotropic), or the measurements represent 
an average over metallographic planes of different random angular orientations and the 
size and orientations are not correlated. In the later case the information concerning 
particle orientations is lost in the averaging process. Metallurgists often come across 
microstructures having particles, whose spatial orientations are anisotropic. The size 
and orientations of these particles effect many of the properties of the materials. Thus it 
is of interest to characterize the ensemble of these particles in three-dimensional 
microstructures to understand their effect on the behavior of materials. In a number of 
cases, there is a strong correlation between the spatial orientation and size of the 
particles. The particles in certain orientations are larger than those in the other 
orientations. For example, in a mechanically worked material, particles that are parallel 
to the compression axis are larger than the particles that are inclined to this direction. In 
such microstructures, the estimation of the joint size and orientation distribution of 
particles provides useful insight into the growth of particles as a function of their size 
and orientation. The sizes and orientations of the observed two-dimensional profiles do 
not necessarily represent the true size and orientations of the corresponding particles 
[72]. Thus, it is of interest to develop a statistical relationship between true bivariate 
particle size and orientation distribution function and the corresponding apparent 
bivariate particle size and orientation distribution function measured on a 
metallographic plane. 

3.3.2 THE MODEL 

In this particular analysis particles are modeled as penny shaped cylindrical 
objects with flat circular ends as shown in Figure 3-11 embedded in a matrix of 
material. When such a system of particles is sectioned with a random plane, the 
intersections of penny shaped particles almost always yield profiles that are chords 
generated by intersections of the flat circular ends with the sectioning plane as shown in 
Figure 3-12. Let the variable R represents the radius of the circular flat ends of the 
particles. Let Rm be the radius of the largest particle. The thickness of any particle is 
expected to be much less than its radius R. It should be noted that this assumption is not 
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necessary for this analysis and the analysis is also applicable to the particles of zero 
thickness. Let the Z-axis of the XYZ coordinate system be the reference axis (or vertical 
axis). Although any direction in the three-dimensional space can be chosen as the 
vertical axis, in practice the direction of applied stress or similar physically meaningful 
direction is chosen as the vertical axis. Any sectioning plane that contains the vertical 
axis is called a vertical sectioning plane (see Figure 3-13). Orientation of a vertical 
sectioning plane is given by the angle (tip(o<<!i^ <2;r). The spatial orientation of the 

particle is given by the angles 0 and (|) pertaining to the vector normal to its flat faces, 
as shown in Figure 3-11, where Q<G<kI2 and 0<(zi^2;T. Let F(R, 9, (^) be the 
joint size and orientation frequency distribution of the particles in three dimensions. 
From the basic principle of geometric probability, F(R, 9, sin9 d9 dcfxiR is equal to 
the fraction of the particles having the size in the range of R to (R+dR), and the 
orientation in the range of 6* to (9+d9) and ^to (<j)+d<f)). Clearly, F(R, 9, (j)) depends on 
the orientation angles 9 and (j). But, it is the orientation of the particle with respect to 
the vertical axis (i.e. angle 9) that is of primary interest rather than the amount of 
rotation around this direction as specified by the angle (f). This is obvious from the way 
vertical axis is chosen. The bivariate distribution F(R, 9) is given as follows. 

Equation 3-63 

2jc 

F(R,9)= \F(R,9,(l>)d(l> 

0 

Consider a vertical sectioning plane of some arbitrarily chosen orientation, say 
<j>p. As mentioned earlier, the intersection of the penny shaped particles with the vertical 
sectioning plane yields traces that are chords as shown in Figure 3-12. Let the variable r 
represent the size (half the length) of these chords. Particles of different radii R yield 
chords of different sizes r ranging from 0 to Rm- The angular orientation of a chord is 
given by the angle cc between the chord and the vertical axis as shown in Figure 3-12. 
Let f(r.a, (pp) be the joint bivariate size and orientation firequency distribution of the 
observed chords on the vertical sectioning plane of orientation in such a way that 
f(r, cc, ^)dr doc is equal to the firaction of the chords having the size in the range r to 
(r+dr) and the orientation in the range cc to (a+da). The fimction f(r, oc, (f>p) can be 
estimated from the experimental measurements of r and a. Let fo(r, a) be the joint 
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bivariate size and orientation frequency distribution of chords, obtained by averaging 
f(r^ cc^ (f>p) over all the vertical plane orientations (i.e. (pp^O to Iji) as shown below, 
Equation 3-64 

1 

fAr,cc) = ~ y[r,a,(l>p)dif>^ 

In actual practice, fQ(r , oc) can be estimated by measuring f(r, cx, ifip) on a few vertical 
sectioning planes of different independent orientations (pp, keeping the vertical axis 
orientation constant and taking an average of these measurements. 

From Equation 3-63, F(R, 6) sinO dR dd is equal to the fraction of particles 
having the size in the range R to (R-^dR) and the angle with the vertical axis in the 
range 9 to (9+d&); the angle (p may be anywhere in the range 0 to 2 tv. Now the 
objective of this analysis is to develop a relation between the true bivariate size and 
orientation frequency distribution of particles F(R, 6) and the experimentally estimated 
apparent size and orientation frequency distribution of the chords fo(r,a). 

The geometry of intersection of a particle with a vertical sectioning plane is 
shown in Figure 3-14. The orientation angle a of a chord in the vertical section of 
orientation (pp, and the orientation angles 9 and (p of the corresponding particle in 3D are 
related as follows, 

Equation 3-65 

sin(^/ip -(p)=CQXa cot 9 

From the geometry shown in (Fig6) the distance between the center of the circular flat 
end of the particle and the vertical plane, d can be expressed in terms of the particle 
radius R, its orientation angle 9, the size of the particle profile observed in the vertical 
plane, r, and the orientation angle of the profile, a, as follows, 

Equation 3-66 

d = 

sma 

It should be noted from the above equation that the distance d is independent of the 
orientation angle (p of the particle and the orientation angle (pp of the vertical plane. 
From Figure 3- 14 it is evident that for a particle of size R r) and the orientation (9, <p), 
as the distance between the particle center and the vertical plane decreases from the 
value d given by Equation 3-66 the profile size increases, but the orientation angle of 
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the profile a doesn’t change as is evident from Equation 3-65. Hence, all the particles 
of size K and orientation (9,(l>), whose center is at a distance less than d from the 
vertical plane, generates chords of size larger than r. Let Ny be the total number of 
particles per unit volume. Let C(r, a, R, 9) be a function such that C(r. a, R, 9, (jff^dR 
d9 da is equal to the average number of profiles per unit area of vertical plane (pp, 
having orientation in the range a to (a-^da) and size larger than r, originating only 
from the particles in the size range R to (R^-dR) having the orientation in the range 9 to 
(9+d9) and (p to ((p+d<p). This quantity is equal to the expected number of particles, 
whose centers lie in a box of unit cross sectional area and thickness 2d {d on either side 
of the vertical plane) as shovra in Figure 3-15. Now we can write. 

Equation 3-67 

C[r,a,R,9,(p^)dR da d9 = ldNy F{R,9,(p)sm9dRd9d<p 

Inserting Equation 3-66 for d in Equation 3-67 we get, 

Equation 3-68 

c{r.a,R,9,(p\dR da d9 = F{R,9,p)sm9dRd9d<p 

^ ^ sma 

From Equation 3-65 it is evident that for given values of ^ 3 nd 9, a and p can’t be 
varied independently. Hence, in Equation 3-68 dp and da are not independent. 
Differentiating Equation 3-65, keeping pp and ^constant we obtain. 

Equation 3-69 

dp cosip. -p) = cot 9 

^ sm a 

Inserting from Equation 3-65 for cos(^p - p) in Equation 3-69 we get. 

Equation 3-70 

dp = cot^ 

Inserting Equation 3-70 for dp in Equation 3-68 yields, 
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Equation 3-71 


c[r,a,R,9,<p^)dR da q a 

sina ^ 


sind dRdOcoiO 


da 1 

^ Vl -cot^ acot^ 9 

^ c[r,a,R,9,cl>^)dR d9 ^ g{a,9)[R^ N, F{R,9,^)dRd9 


Where 

Equation 3-72 


g{a,e) = 

sin’ Q:[sin^ a -cos^ dY^ 

Let Co(f, cx, R, 9) he the value of C(r, a, R, 9, (ftp) averaged over all the vertical plane 
orientations, keeping r. a, R, <9 constant. Then we can write, 

Equation 3-73 

*1 2/r 

C,(r,a,R,d) = — jc{r,aM</>p)d(pp 

IK Q 

To calculate this average, integrating Equation 3-71 over gives. 

Equation 3-74 

]c(r,a,R,9J^)d<p^dR d9 = g{a,9)[R^ -r^f^ Ny ]F(R,9,(t))d(l)/Rde 
0 0 

In the above integration, as mentioned earlier r, a R, 9 are to be kept constant. Thus, 

only the orientation angle ^ of the particle may vary as the vertical plane orientation ^ 

changes. It should be noted at this point that the expressed by Equation 3-70 is at 

constant <^p and 9, whereas the d(f) we are considering at this context is at constant a, 9. 

At constant a, 9, Equation 3-65 dictates that d<l> must be precisely equal to d(f>p. Thus, 

we can intepute F(R, 9, (j>) in Equation 3-74 over instead of Now, Equation 3-74 

can be rewritten as 

Equation 3-75 

=s(a,B)[R^ -r'f N, 

0 ° 

Inserting Equation 3-63 and Equation 3-73 in Equation 3-75 yields. 
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Equation 3-79 




ArN^ 


/r /2 


/r A/’^ sin^ a 


J J 

rixll 


cos- 9^\ae F{R,e)dedR 
sin^ a — cos^ 9 


The above double integral equation can be changed into the form of standard double 
Abelian integral equation and can be solved in the same manner as described in section 


3.2.2. 1 to give the following general solution. 
Equation 3-80 


^ ffjl p)- "f'V sin^ acosaf,{r,a)dadr 

^ ^ ' 4cos^^sin^]5i25^ J j [cos^- a - sin^ ofir^ ~ R^P 

The quantity on the left hand side of the above equation is the three dimensional 
counter part of nA(r, a) and is denoted by ny(R, 9). Thus the quantity ny(R, 9) sin9dR 
dO is equal to the number of particles per unit volume of the three-dimensional 
microstructure having size in the range R to (R+dR) and the orientation in the range 9 
to (&+dd). Now the above equation can be rewritten to give the final form of the 
expression for estimation of the tme bivariate particle size and orientation distribution 
from the corresponding measured average apparent bivariate distribution on the vertical 
sections as. 

Equation 3-81 

(d a'\ - 1 sin^g cos g/t^(r,ayQair 

’ 4cos^^sin6»]5i?5^ j j [cos'a-sin'^}^[r' 

It should be noted that none of the above equations contains the particle thickness term 
and hence all of the above equations are applicable for any thickness value of the 


particles, including zero. 








Figure 3-5 Polydispersed system of spheres 
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Figure 3-7 Schematic illustration of the sectioning of spheres of diameters Di to D5 to 
give circular section of diameters 2ri to 2r5 
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Figure 3-8 Transverse and longitudinal section of spheroid. Axial ratio of spheroid is 
reflected in axial ratio of longitudinal section. 
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PROUTE SPHEROIDS 



/ 

OBLATE SPHEROIDS 


Figure 3-^) Shape variation by variation of axial ratio within the types of spheroids 
(prolate or oblate) 
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Figure 3-1 1 Model shape for microcracks, and specification of microcrack onentation. 
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Figure 3- 1 2 Specification of size and orientation of microcrack traces in a vertical 

plane. Note that the vertical axis is contained in all the vertical planes, and 
hence it can be identified in each vertical plane. 


73 



Figure 3-13 


(a) 


f Vertical axis 




(a)Vcrtical planes of different angular orientations (b) Vertical planes at 
different locations in the specimen and with different orientations. 






c 



()(.)' i lie 

OQ 1 vertical plane 

O' ; center of microcrack trace 

I'igurc 3- 1 4 ( icomctry involved in the intersection of a penny shaped microcrack and a 
vertical plane. A decrease in the distance d brings the microcrack closer to 
the vertical plane and increases the size r of the microcrack trace. 


75 




CHAPTER 4 


RESULTS AND DISCUSSION 

4.1 ANALYSIS OF NUMERICAL METHODS 

4.1.1 COMPARISON WITH EXACT SOLUTIONS 

The numerical procedures developed are validated for their exactness by using 
raw data, for which exact solutions are known before hand. Basically computer 
programs are written to solve the problem of numerical de-convolution of 2D data to 
obtain the 3D distribution for the following three categories of particles; 

a) Monodisperscd tmd polydispersed systems of spheres 

b) Prolate and oblate spheroids of varying size and shape 

c) Platelets having varying size and orientation. 

In each case raw ilata is developed from a known 2D distribution function, for which 
exact 31) solution is known. 

4. 1.1.1 SPHERES 

In case of spherical particles, the 3D probability density flmction ,F( R ) is 
assumed to be a nonnal distribution function with a mean value of 10 and a standard 
deviation of 3. The corresponding probability density function for expected circular 
profile sizes ll(r) is obtained by solving this function analytically. The probability 
densities of both of these functions are shown in Figure 4-1. A universe of 5000 
profiles is generated using f(r). 

The generated 2D data is then fed into the computer program developed to 
estimate back the 3D particle size distribution and 3D mean particle size, and the 
results obtained are presented in and Figure 4-3. 3D particle size distributions are 
obtaincil using H, 1 5 and 25 number of size classes. As can be seen firom, although the 
sample size is large enough (5000) the estimated 3D distribution is not perfectly 
matching with the exact 3D distribution when the particle sizes are descretised into 8 
classes. The error is not because of small sample size (since 5000 is a large enough 
sample size) but because of using inappropriate descretisation of particle sizes. This is 
evident from figures showing estimated 3D size distribution using 15 and 25 size 
cla.s.ses ( {b),(c)), where the estimated 3D size distribution matches well with the exact 

3D size distribution. 
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From the estimated 3D particle size distribution, mean 3D particle size is calculated in 
all the three cases and also in cases of 10 and 20 size classes and the results obtained 
are presented in figure 4-3. It is evident that the estimated mean size approaches the 
exact value on increasing the number of bins. This is because on increasing the number 
of bins the class interval decreases providing more detailed description of the 
distribution, Altlu)ugh this is only applicable if we have a large sample size. It should 
also be noted that the mean profile size (2D) is less than the real 3D mean particle size 
in this particular case. In general, the mean profile size can be either less, more or equal 
to the real 3D mean particle size (e.g. bimodal distribution of 3D particle sizes). Hence 
using the mean profile size as a representative of real mean 3D particle size may lead to 

erroneous results, 


4.1.12 VARlABIJiSinUvROlDS 

Applicaliility of the numerical solution developed for the estimation of 3D 
bivariate si/c-shape distribution of variable spheroids is tested. This involved de- 
convolution of 2D bivariate size-shape distribution of corresponding elliptical profiles, 
generated on 2D sectioning plane by spheroids. A hypothetical population of variable 
prolate spheroids, whose size and shape are governed by the following joint probability 
density function, (71 j is considered. 


Equation 4-1 

l’{h,x’)==( - *“ )+ (l - ^’)1 

The three-dimensional representation of the above equation can be seen in Figure 
4-4(a). The joint probability density function of the coiresponding elliptical profiles is 

obtained Iroin lirioation ,l-14andisgivenby 

liquation 4-2 

-/Kl-wO+(l-/)] 




■I bis is graphically represenled In Figure 4-4(b). Using this 2D bi«r«fe s.e-^^e 
probahilily densily funolion a universe of 5000 ellipfical profiles ts genera^d. e 

siae of minor principal axis varying from 0,01 fo 1.0 and dte square ec^ntnety 

parameter varying from 0 to I . 

The gcneralcd 21) elliptical profile data is then fed uito t = 

devel„p«lloes.i,..nte.he51)blvariate.z.shapedisulhutlonofv.^^^^^^^ 

the results obtained are presented in Figsrre 4-5 and F.^rre 4-6. Two case 
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considered for this analysis. In the first case the particle sizes and shapes are 
descreti.sed into agiid ol 10 X 10 size. Whereas in the second case a size-shape grid of 
25 X 25 is used. Similar to the sphere size distribution problem, the effect of class 
interval on the. estimated 3D distiibution can be seen froni Figure 4-5(a-b) In the first 
case the estimated 3D bivariate size-shape distribution of variable prolate spheroids is 
significantly dilTei ent from the exact 3D bivariate distribution, although the sample size 
is quite large {5000). Whereas, in the second case the estimated and exact 3D bivariate 
size-shape distiibuticms arc almost similar. The histogram representation of the exact 
and the estimated 3D size-shape distributions are represented in Figure 4-6. 

4.1.1 .3 PI A'W.l .in s 11 AVINCi VARYING SIZE AND ORIENTATION 

A computer simulated population is generated to validate the suitability of the 
numerical piocedme developed for the estimation of 3D size and orientation 
distribution of disk shaped particles (platelets). The population consists of rectangular 
traces genciateti from hy[Hnhctical joint probability density of size and orientation 
governed by the (bllowing expression (72] 

Equation 4- 1 

[l2«10’O4|r.)“-(ry] 

The underlying assumption in solving the size orientation problem leads to a fact that a 
randomly oriented (in this context) disk is not actually randomly oriented with respect 
to the solid angle distribution The joint probability density function for corresponding 
disks in 3D can be obtained from Equation 3-19and is given by 
Equation 4-2 

/r, {R,0) 64x 10“/|/eJ^ 

Both the above exact joint probability density functions (Equation 4-1 for 3D disks and 
Equation 4-2 lor concsptmding 2D profiles) for 3D discs and are graphically 
represented in Figure 4-7 A universe of 5000 2D rectangular profiles is generated 
using Equation 4-1, with the size of the traces varying from 0.00005 to 0.001 and the 
(mentation from 0 to 00 degrees, 

Estimation of 31) bivariate size-orientation distribution of disk shaped particles 
is generated by nutnci ical de-convolution of the above 2D data. The results obtained 
arc prcsentetl iti Figure 4-K and Figure 4-9. The size and orientations of the particles are 
descretised into 25 classes each. It can be clearly seen from Figure 4-8 and Figure 4-9 

that the estimated tr ue bivariate size and orientation distribution of disks almost follows 

79 



,he corresponding exact bivariate distribution, except for a negligible discrepancy in the 

frequency values. 

4.1 .2 ANALYSIS OF APPROXIMATION METHODS 

The objective of this analysis is to study the appropriateness of assuming 
spherical shapes for non-spherical particles to de-convolute 3D particle size 
distribution. For this analysis (he same universe of five thousand elliptical profiles as in 

4. 1 . 1 . 2 is consiilered. I he numerical procedure developed for 2D to 3D conversion of 
spherical particles is applied to estimate the size distribution and mean particle size of 
this non-spherical system. This analysis basically consists of two parts. 

a) Arithmetic approximation 

b) (lemnetric approximation. 

In arithmetic apiproximation, an elliptical profile is approximated as a circle of radius 
‘r’ equal to the arithmetic mean of the principal semi axes of the elliptical profile. Thus 
an elliptical profile of minor principal semi axis ‘m’ and square eccentricity parameter 
'y' is approximated a.s a circle of radius T’ where, 

Equation 4-5 



In case of geometric approximation, an elliptical profile is approximated as a circle of 
equivalent area. I'lnis an elliptical profile of minor principal semi axis ‘m’ and square 
eccentricity parameter ‘y’ is approximated as a circle of radius T’ where, 

Equation 4-6 



ILsscntially luiualion 4-5 mul Equation 4-6 assume that the radius of circular profile is 
an arithmetic mean ami a geometric mean of the two elliptical radii respectively. 

The gcnci atctl universe of 2D elliptical profiles is converted into 1) universe of 
equivalent circular profiles u.sing arithmetic approximation 2) universe of equivalent 
area circular profilc.s u.sing geometric approximation. Then the corresponding 3D size 
distribution of particles is estimated for both these universes by using the numerical 
routine deployed in 4. l.I.l. Since the estimates provide size distributions of spheres, 
whereas the tnie 3D di.siribution is for bivariate spheroids, we need to convert the true 
3D distribution to an equivalent spherical distribution. This is again done using 
arithmetic and gct)meiric approximation by taking equivalent sphere radius as 



arithmetic mean and geometric mean of radii of spheroids. The results obtained are 
presented m lagure 4-10, Figure 4-11 and Figure 4-12. The error involved in 
approximating the non-spherical particles as spherical particles in estimating the trae 
particle .si/.c di.sinhution is evident from Figure 4-10. It is observed that the geometric 


approximation g.ivais better results than arithmetic approximation. In Figure 4-10 we 
can clearly ace that the cui-ve representing the estimated true particle size distribution in 
case of I'comctnc approximation follows almost the same trend as the curve 
representing the exact true particle size distribution, except for some significant 
discrepancy in the frequencies, whereas the curve representing the estimated true 


particle .size ilistribution in case of arithmetic -approximation is more off from the exact 
true particle .si/e distribution. However, in both the cases, it is observed that the 
approximation mclhod always results in overestimation of the number of small particles 
and underestimation of large particles. sTiows the same graphs shown in Figure 4-10 but 
arc wiiicly spread near the large particle size range. It can be seen in Figure 4-1 1 that 


the approximation of non-spherical particles as spheres also leads to prediction of some 
large fiarticlcs, which arc actually not present in the true population of variable 
splicroid.s, 

The mean true particle size of the non-spherical system under consideration is 
also cstimatctl from the .size distribution obtained by approximation methods. The 
results ivhtained arc depicted in Figure 4-12, In this figure three average sizes are shown 
for each ca.se of approximation. The exact mean size corresponds to the mean of 
equivalent .spherical radius of the true population of variable spheroids. The bivariate 
mean size is similar to the exact mean size except that the bivariate spheroid 
disirihution i.s not Inic distribution but has been generated by de-convoluting 2D 
bivariate distulmtion as outlined in 4. 1.1.2. The approximate mean size corresponds to 
the mean size ol' the distribution of spheres which is generated from 2D ellipses 
a.sKummg them a.s circles. It can be seen from Figure 4-12 that the geometrical 
approximation method gives better result than the arithmetic approximation method, fri 
fact in thi.s (larticular ca.se the geometrical approximation method gives almost the same 
result a.s obtained by the bivariate estimation method. It is interesting to note that 
although the approximation methods overestimate the smaller size particles and 
undcrc.stimutc the large size particles, they give a mean particle size that is more than 
the exact mean particle size. It may be because of the fact that the approximation 

mcthod.s estimate some non-exi.sting very large size particles. 
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4 #2 SAMPLING REQUIREMENTS AND ASSOCIATED ERRORS 

The objective ol- this analysis is to emphasize the importance of sample size in 
the estimation of 3D particle size distribution and 3D mean particle size by numerical 
dc-, convolution of apparent 2D profile size distribution. For this analysis the same 3D 
probalnlitv distribution function of spheres deployed in 4.1. 1.1 is used to generate a 
universe of lOOOO. 2D circular profiles. From this universe 5 samples each of sizes 100, 
200, 500, 1000 and 2000 are taken at random. These samples are used to generate 
estimations of 3D di.stribution and are stored as a data set. 

The numerical procedure developed for de-convolution of 2D data to 
approximate the true sphere size distribution is run with all these randomly generated 
datasets. In all thc.se cases the particle sizes are descretised into 8 classes. The results 
obtained arc presented in I'igurc 4-15. It is seen from Figure 4-15 that the estimated 
piirticlc si/c distribution approaches the exact particle size distribution with increase in 
the sample si/.c. Tlii.s is due to decrease in sampling error with increasing sample size. 
It may be noted that this seems similar to the predictions of central limit 
theorem! uiuicr wood I, though the theorem is not exactly applicable since this not a 
sampled raiulom variate. It is also noticeable from Figure 4-15 that the precision 
obtained in the estimation of particle size distribution is not a linear function of the 
sample .size, The cstimatctl piirticlc size distributions are significantly different from the 
exact piirticlc si/.c di.stnhution in case of samples with sizes of 100 and 200. Whereas, 
samples of si/.cs 500, 1000 and 2000 are all very close to the exact distribution with 
only minor improvement.s in estimation with increasing sample size. This implies that 
increase in .ssimplc si/.c don’t yield corresponding betterment of estimates beyond a 
certain sample size. 

The effect of sample size on the estimation of mean 3D particle size is also 
analyzed and the lesults obtained arc presented in Figure 4-13 and Figure 4-14. The 
estimatcii mean values show a lot of scatter around the true mean, which decreases with 
increase in sample size. Tliis decrease in scatter with increasing sample size is obvious, 
.since it dccri.‘a.scs sampling error as stated in previous paragraph. Further, it is noted 
that the scatter is not uniformly distributed about the true mean and is predominantly an 
over estimate rather than random, for all the sample sizes. This is a very peculiar 
observation, imliculing a bias in the process of de-convolution, since the sampling was 
randomly done without any bias. It was realized that the mathematics of de-convolution 
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w<is cxuct <in(i thus tfic- tinly source ot bias could be due to the numerical method. This 
meant that either the bi.is vviis due to descretisation or floating point operations. Since 
the bia.s wa.s appreciable it was concluded that it couldn’t be due to 8-bit floating point 
operation.s. 1 Icnce by process ot elimination it was assumed that descretisation was the 
most probable reason toi the bias. On revisiting the numerical method for descretisation 
it was reali/ctl tliat the size class is being indicated by the arithmetic mean of the 
interval. .Since the bias was ;m over estimation and it is a well known fact that 
arithmetic mean is greater than geometric mean and harmonic mean, it was anticipated 
that representation of mean bin size by geometric mean and harmonic mean would 
result in lowering the mean probably removing/decreasing the bias. The results of 
application of i'ct^tH'bic mean and harmonic mean are shown in Figure 4-14. As 
cxpcctcii the estimates of mean size are lower than that of arithmetic mean. Although 
these seem to he al.-.o hiaseil with an under estimate. Hence it was realized that a more 
appropriate fniietion shouUl be used to represent the size class in the process of 
descretisation. As a tc.st function root mean square value of the interval was taken as 
mean si/.e ot the bm. Ihc calculations based on this are shown in Figure 4-14. It is 
interesting note tlnit use of root mean square value for mean bin size results in an 
evenly disti ibulca .scatter around the exact mean, indicating absence of any bias. It is 
also seen tliat the .spread in the scatter si^ificantly decreases on using the root mean 
square vtilue. I icncc, root mean square value should be the function of choice m 
representing the mean bin si/.c in solving the de-convolution problem for spheres. 



Exact Distribution Densities 



Figure 4- 1 A splierc size tlistributioii and the corresponding profile size distribution. 

Solid line corre.sponds to the sphere size distribution density , whereas the 
dotted line represents the probability density of corresponding circular 
profile .si/cs. 
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Figure 4-2 Comparison of the estimated 3D sphere size distribution with the exact 3D 

sphere size distribution. The apparent 2D circular profile size distribution 
is also shown. In a) Particle sizes are descretised into 8 classes, in b) 
Particle sizes are descretised into 15 classes and in c) Particle sizes are 
descretised into 25 classes. Same 2D profile data is used in all the three 
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cases. 
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Figure 4-3 Illustration of variation of estimated 3D mean size of the spheres with 
change in number of size classes. Notice that the estimated mean size 
moves closer to the exact 3D mean size with decreasing class interval. • 



Exact Distribution 



2D Size Shape Distribution 



iMgurc 4-4 3D representation of size-shape probability density function. The surface 
in (a) describes a hypothetical population of variable, prolate spheroids. 
The surface in (b) corresponds to the variable elliptical profiles 
determined in the spheroids by random sectioning. 
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Exact vs Estimated Distn(25 Bin) 



Exact vs Estimated Distn(25 Bin) 
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[■'igurc 4-5 Comparison of the estimated 3D size-shape probability density of variable 
spheroids (represented by Figure 4-4(a)) obtained by numerical de- 
convolution from 2D data with the exact 3D size-shape probability 
density. In (a) Particles are distributed into 10 size classes and 10 shape 
classes and in (b) Particles are distributed into 25 size classes and 25 shape 
classes. 
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Exact 3D Distribution 



4..(, Histogram representation of 3D size-shape distribution of (a) Population 
of variable spheroids shown in Figure 4-3(a), (b) Population estimated by 
numerical de-convolution with a 10 by 10, size-shape grid and (c) 
Hstimated population with a 25 by 25 size-shape grid. 
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Figure 4-7 3 D representation of size-orientation distribution function. The surface (a) 

describes a hypothetical population of randomly oriented disks. The 
surface in (b) corresponds to the apparent disk trace distribution in vertical 

planes. 
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Exact vs Estimated Distn(25 Bin) 



igure 4-8 Comparison of calculated tme bivariate size and orientation distribution 
obtained by numerical de-convolution of apparent bivariate size and 
orientation distribution of microcrack traces (represented by Figure 4- 
6(b)) with the exact size and orientation distribution of microcracks. 
Particle sizes are descretised into 25 classes and orientations too are 
dcscretiscd into 25 classes. 
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Exact Vs Estimated Distn. with Geometric 
Approximation 



I‘igurc4-1() Spurious distribution curves resulting from application of spherical 
conversion coefficients to non-spherical particles (prolate spheroids, in 
this case). Solid lines are the exact 3D size distributions. In (a) the 
elliptical profiles are approximated as circles of radius equal to the 
arillunelic mean of their major and minor principal semi axes. In (b) the 
elliptical profiles are approximated as circles of equal area, i.e circles of 
radius equal to the geometric mean of their major and minor principal 


.seini axes. 
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'I'hc same distribution curves as shown in Figure 4-10 are exploded to 
show tlic difference between the estimated and exact distribution curves in 
the larger particle size region. 
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I'ij,>.ur(.‘ 4-12 Comparison of mean particle size obtained by approximating the non- 
sphcrical particles as spheres (both arithmetic and geometric 
approximation) with the exact mean particle size. Mean particle size 
olMaincd by using the actual non-spherical conversion coefficients is also 
shown.' 






Figure 4»13 Illustration of variation of estimated mean particle size with the sample 
size (100, 200, 500, 1000, 2000). In (a) size of the particles in a particular 
bin is taken as the arithmetic mean of the sizes corresponding to the ends 
of the bin. In (b) Geometric mean is taken as the particle size. In (c) 
Harmonic mean is considered as the particle size and in (d) root mean 
square value is taken as the particle size. 








Figure 4-14 Illustration of bias involved in considering the particle size as ari 

mean size, geometric mean size, hamonie mean size and toot mean 
square size in estimating the 3D mean particle size by numerical de- 

convolution of 2D data. 
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4- 1 5 Eiffect of sample size on the precision, 

estimation of 3D particle size distribution by numerical de-convolution of 
2D data. 5 samples each of sizes 100, 200, 500, 1000, 2000 are considered 
lor this analysis. 








CHAPTER 5 


CONCLUSIONS 

This study ot 2D to 3D de-convolution problem led to the following 

conclusions. 

1 . Numerical methods are developed for de-convolution of observed 2D profiles of 
spheres, variable size-shape spheroids and variable size-orientation disks were 
successfully developed. 

2. The numerical methods for solving the do-convolution problem for all the three 
(spheres, variable spheroids, variable disks) yield very precise 3D distributions 
(provitied enough sampling is done). 

3. For large sample size fine descretisation is required to achieve better estimations. 

4. Use of profiles of spheroids as equivalent area circles to de-convolute 3D size 
distribution is not a bad approximation. 

5. Approximation of spheroids as spheres led to prediction of unrealistic large 
particles and also increase frequency of lower size class particles. Though the two 
factors negated each other in the estimation of mean. 

6. The error in estimation can be decreased by increasing sample size. Though the 
increase in sample size don’t yield corresponding betterment of estimates beyond a 
certain sample size. This sample size was seems to be about 500. 

7. In de-convolution of spheres the descretisation process laeads to a bias in the 
estimate. It is a positive bias when the mean bin size is represented by arithmetic 
mean of interval and a negative bias when it is represented by geometric mean and 
harmonic mean of the interval. 

8. Use of rot)t mem square value (in de-convolution of spheres) to represent mean bin 
size provides the most bias free estimate with least spread about the exact mean. 
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CHAPTER 6 


(ailDELINES FOR FUTURE WORK 

1 ) The numerical methods developed in this thesis can be further analyzed for their 
applicability with other kinds of distributions like lognormal distribution, 
bimodal distribution etc. 

2) Effect of standard deviation of the exact distribution on the accuracy of 
estimated distribution can be studied. 

3) Numerical method for unfolding trivariate size, shape and orientation 
distribution of spheroids will be a more helpful microstructure characterization 
tool, since spheroidal particles often have a trivariate distribution. 
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